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ABSTRACT 

Data from the Hubble Space Telescope taken with the ACS/WFC camera have been used to in- 
vestigate the globular cluster (GC) populations around six giant elliptical galaxies that are ~ 40 
Mpc distant. From these six fields, imaged in B and /, a total of more than 15000 candidate GCs 
have been measured, of which 8000 or more are high-probability globular clusters. The data reach a 
limiting magnitude near Mj ~ —8, about 0.4 mag fainter than the GC luminosity function turnover 
point, and thus thoroughly cover the bright half of the GC population. Most of the individual GCs 
on these images are marginally resolved nonstellar objects, so King-model profiles convolved with the 
stellar point-spread functions are used to measure their individual total magnitudes, colors, and linear 
effective radii. The classic bimodal form of the GC color- magnitude distribution shows up unam- 
biguously in all the galaxies, allowing an accurate definition of the mean colors along each of the two 
sequences as a function of magnitude (the mass/metallicity relation or MMR). The blue, mctal-poor 
cluster sequence shows a clearly defined but nonlinear MMR: in this particular photometric dataset 
the mean GC color changes smoothly from a near-vertical sequence at low luminosity {Mj > —9.5) to 
an increasingly redward slope at higher luminosity. By contrast, the red, metal-rich sequence shows 
little trace of a MMR and is nearly vertical at all luminosities. The form and slope of the MMR along 
either sequence do not depend strongly on either cluster size or galactocentric distance Rgc- 

All the observed features of the present data agree with the interpretation t hat the MMR is create d 
primarily by GC self- enrichment, along the lines of the quantitative model of iBailin fc Harri^ ()2009l) . 
During the protocluster formation stage, the more massive GCs are better able to hold back the 
enriched products of the earliest supernovae and to seed the lower-mass stars still in formation. The 
"threshold" mass at which this effect should become noticeable is near 1 million Solar masses, which 
is closely consistent with the transition region that is seen in the data. More generally, the data favor 
models in which the star formation efhciency in a protocluster is roughly independent of mass, and in 
which the gas retention efhciency is a strong function of mass. 

Correlation of the median scale sizes r^ of the GCs with other parameters shows that the metal- 
poor clusters are consistently 17% larger than those of the metal-rich clusters, and that this difference 
holds at all galactocentric distances and luminosities. At the same time, cluster size scales with halo 
location as r^ ~ ^^'c^^ indicating that both metallicity and the external tidal environment play roles 
in determining the scale size of a given cluster. Lastly, both the red and blue GC components show 
metallicity gradients with galactocentric distance that are shallow but real: heavy-element abundance 
scales as Z ^ ^gc*'^ fo'" both types. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: star clusters - globular clusters: 
general 



1. INTRODUCTION 

Much of the recent observational work on the pop- 
ulations of globular clusters (GCs) in other galaxies 
has concentrated on their metallicity distribution func- 
tion (MDF). A near-universal result, confirmed from 
both photometric and spectroscopic samples of GCs in 
many galaxies of all types, is that the MDF is bimodal, 
strongly suggestive of two major epochs of cluster forma- 
tion (see, e.g. Zcpf & Ashman 1993: Larson ct al. 200i!; 
Harris et a l. 2006; Peng ct al. 2006; Kundu & Zcpf 2007; 
Strader et al 2007 : Wehner et al., 2008 : Pen g ct al. ,200a; 
Waters et al.ll2009l for an extensive range of results and 
discussions). This bimodality paradigm has steadily 
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been reinforced as the quality and size of the databases 
have increased. 

A more recent discovery is the existence of intriguing 
second-order structure in the MDF. In some large galax- 
ies a correlation between color and luminosity has been 
found along the blu er, more meta l-poor G C sequence 
(iHarris et all l2006t iStrader et all [2b06t iMieske et all 
I2006i) . Since the integrated color of old star clusters de- 
pends predominantly on their metallicity, this correla- 
tion represents a mass/metallicity relation or MMR (the 
trend is also referred to colloquially as a blue tilt, but 
"MMR" is more accurate and more general) . The essence 
of the trend is that the more massive blue GCs are pro- 
gressively more heavy-element enriched. Just as intrigu- 
ingly, the redder and more metal-rich GC sequence has 
not shown any definite evidence for a MMR even in galax- 
ies where the blue-sequence MMR is strongly present. 
Early hints o f this e ffect were, in fact, clearl y noticed by 
lOstrov et al! (|1998f ) and lDirsch et al.l (|2003f l from wide- 
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field photometric surveys of the GC system in the Fornax 
cD NGC 1399. There, they pointed out that the highest- 
luminosity GCs showed a broad unimodal color distri- 
bution, as if the blue and red components had merged 
together at the top end. 

To make the effect more challenging to understand, 
some large galaxies su ch as NGC 4472 may not show the 
blue- sequence MMR (iStrader et al.l [200a [MI eske et al.l 
[20061) even though the measured sample of GCs seems 
large enough to re veal it, if present. M 87, at first claimed 
to have a MMR (jStrader et al.ll2OO60 . doe s not show it 
accord ing to a larger set of data analyzed by lWaters et al.l 
(|2009D . As will be discussed more extensively below, a 
critical factor in deciding the presence of an MMR is not 
just the total measured population of GCs, but also the 
highest luminosity (mass) to which the sample reaches. 
These factors are also crucial in deciding whether or not 
the MMR might be linear in color versus magnitude (as 
is assumed in most previous papers) , or nonlinear, which 
in turn is a central issue for the model interpretations 
that are now starting to appear. 

Early descriptive interpretations for the physical cause 
of th e MMR have been suggested in. the discovery pa- 
pers ([Harris et al.ll2006l : IStrader et a"ni2006t iMieske et all 
|2006[ ). but the most promising direction at present for 
producing the e ffect is organized around the idea of GC 
self-e nrichment (jStrader fc SmithI 120081 : iBaihn fc Harris! 
|2009( ). The basic approach is that clusters forming within 
massive protocluster gas clouds can hold on to some por- 
tion of their SNe ejecta during the first round of star 
formation and thus enrich the still-forming lower-mass 
stars in the cluster. The higher the proto-GC mass, the 
higher the gas retent ion, thus leading to a MMR. The 
quantitative model of iBailin fc Harrig (|2009( l shows that 
the amount of self-enrichment should be negligible for 
protoclusters with M < lO^Af©, which are small enough 
that SN ejecta easily escape their potential well. Thus for 
clusters smaller than this threshold, the cluster metallic- 
ity should be uncorrelated with mass; that is, the GC 
sequences should be vertical below this transition-mass 
point. 

In addition, the Bailin/Harris model predicts that a 
red-sequence MMR should exist, though this should be- 
come noticeable only at still higher GC mass and with 
a smaller amplitude because the effects of internal self- 
enrichment will be a smaller fraction of the initial (pre- 
enriched) cluster metallicity at any mass. There is even 
a hint that this red-sequence MMR has already been ob- 
served in the extremely populous cluster systems in NGC 
3311 and NGC 4874 (|Wehner"eiran [20?)8l: iHarris et all 
[20091 : IBailin fc Harrisll20(W ! 

At the same time as the interpretive models have ad- 
vanced, the state of the obser vations is stil l some what 
confused. The initial claim of [Harris et al.l ([2006( ) was 
that, at masses less than ~ IO^Mq, the blue sequence is 
vertical (i.e., color is uncorrelated with luminosity) and 
that only the higher-mass range showed the correlation. 
Other papers discuss the MMR in terms of a linear re- 
lation between c olor and metallicity extending continu- 
ously downward fStrader et al.' 2006; 'Spitler et al .|[2006[: 
CanticUo ct al. 2007; DcGraaff ct al. 2007; Wchnci{et3] 
20081 : ISpitler et al.ll2008t iCockcroft et al.ll2009l) . 

Another, and extremely important, source of confu- 
sion is understanding why various individual galaxies do 



not show an obvious MMR; for example, comments have 
been raised that the Milky Way does not have the effect. 
However, the Milky Way does not have enough clusters 
to reveal a MMR with any statistical confidence whether 
or not it is present. The effect of small-number statis- 
tics is especially important if the MMR affects only the 
higher-mass GCs, which are the rarest. 

With the benefit of hindsight, it is easy to under- 
stand why these correlations were not noticed earlier: the 
MMR slope is relatively modest, so big statistical sam- 
ples of clusters and precise, metallicity-sensitive photom- 
etry are both needed to see it. Few such samples existed 
till quite recent years. One way to put this into perspec- 
tive is to note that the basic phenomenon of bimodality 
is a first- order feature of the GC metallicity distribu- 
tion that becomes obvious even with samples as small 
as ^ 10^ clusters (the prime example is the Milky Way 
itself, which was the first galaxy in w hich the bimodality 
paradigm became clear (|Zinnl[l985[ )). By contrast, the 
MMR is a second-order feature and thus in rough terms 
can be expected to require samples of ^ 10'^ clusters to 
become clearly detectable. 

In giant ellipticals it is possible to measure hundreds 
of GCs well above the million-Solar-mass level where the 
MMR is agreed to be most visible; for example, the GC 
sample to be described in this paper (see below) has 750 
clusters brighter than Mj = —10.5, whereas in the Milky 
Way there is only one such cluster, oj Centauri. The 
same comparison is true of dwarf galaxies and most disk 
galaxies, where any attempts to search for the MMR are 
severely hampered by small-number statistics even with 
high-quality photometry or spectroscopy. In other words, 
for most galaxies on an in dividual basis, the ex istence of 
a MMR is not decidable. [Mieske et all ([20061 ). in their 
study of GCs in 76 Virgo galaxies, worked around this 
issue by combining their galaxies into four luminosity 
groups (dwarfs through supergiants) to gain increased 
statistical weight. Even so, the dwarf sample contains 
too few GCs in the critical high-luminosity regime to 
clearly reveal a color/luminosity slope there. 

A challeng e to the e xisten ce of the effect itself has been 
presented bv lKundd ([20081 ). His claim is that the MMR 
is spurious, an artifact of the combination of aperture 
photometry, low signal-to- noise, and a mass/radius re- 
lation for GCs. These concerns will be dealt with in 
this paper in the sections below. It is, however, al- 
read y worth em phasizi ng that the three discovery pa- 
pers ([Harris et al 2006. : .Strader et al.l[2006l : [Mieske et all 

[2006f ) employed three different photometric measurement 
procedures over a wide variety of galaxies, yielding simi- 
lar MMR slopes. In addition, ground-based photomet- 
ric studies of the GCSs in distant galaxies, in which 
the individual GCs are completely unresolved and any 
differential-size issues become moot, also show the MMR 
([Dirsch et al.[ [20031 : [Wehner et all [20081 ) . Nevertheless, 
given particularly the growing state of the theoretical 
modelling, the data need to be put onto a better foun- 
dation. 

The largest GCS populations of all, and in principle the 
best for a clear measurement of the MMR, are found in 
giant and supergiant ellipticals, many of which are cD's 
or Brightest Cluster Galaxies (BCGs) at the centers of 
rich c lusters ([Harris! [20011: [Harris et all [1995!: iBlakeslei 
[1991 IHarris et al.ll2006[ ). In this paper, I present new 
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and homogeneous photometry for six gE's with globular 
cluster systems that are among the largest known. 

Throughout this paper, a distance scale Hq = 70 km 
Mpc~^ is used to convert redshifts and angular di- 
ameters to true distances and linear radii. 

The plan for this rather long paper is as follows: Sec- 
tion 2 introduces the set of galaxies being studied and the 
goals for the measurement, while Section 3 steps through 
the complete process of measurement of the cluster mag- 
nitudes, colors, and linear sizes, and Section 4 presents 
an overview of the final color-magnitude data. Section 5 
describes the complete analysis of the GC effective radii 
(scale sizes) and traces their correlations with external 
parameters such as metallicity, galactocentric distance, 
and metallicity. Section 6 presents the full results for the 
mass/metallicity relation along both the blue and red GC 
sequences, and its implications for the theoretical mod- 
els. Lastly, Section 7 presents new results for the spatial 
metallicity gradients in the red and blue GC subsystems. 
The discussion concludes with a summary of the findings 
in Section 8, and some concerns in Section 9 that may 
affect both the basic nature of the observations and its 
links with theory. 

2. THE DATA SAMPLE 

The raw data are for six galaxies drawn from the giant 
ellipticals observed in HST program GO 9427 (Harris, 
PI) . These six, listed in Table 1 , are all imaged with the 
ACS/WFC in the two broadband filters FA35W (B) and 
F814W^ (/), and are the ones in this particular program 
that turned out to have the largest GC populations. As 
noted in Table 1, most are at distances ~ 40 Mpc. The 
data from this program also contained three other gi- 
ant E targets, NGC 5322, 5557, and 7049: these three 
had distinctly smaller GC populations, and as will be de- 
scribed below, they were used instead to help define the 
level of field contamination.^ 

In the origina l data paper for this program 
(jHarris et al.ll2006l ). both PSF- fitting and fixed-aperture 
photometry were used to measure the GCs. However, 
neither approach is strictly correct, because many of the 
GCs are not quite small enough in linear size to appear 
starlike on the images, and furthermore their intrinsic 
sizes are not identical from cluster to cluster. Various 
tests in the literature (e.g. iLarsenl 119991: iJordan et al.l 
[20051: iMieske et al.ll2006l: iHarris et al.ll2009l) demonstrate 
that when the FWHM of the star cluster is smaller than 
about 10% of the FWHM of the stellar point-spread 
function (PSF), its size can no longer be determined 
from these conventional imaging techniques and it can 
be treated as starlike. But if it is bigger, its size needs 
to be individually accounted for to determine a correct 
total magnitude. For purposes of comparative discus- 
sion, we can usefully distinguish four different regimes of 
resolution: 

• Well resolved: The FWHM or effective diam- 
eter of the GC is much larger than the PSF, 

^ NGC 7626 was not included in the 2006 paper but is added 
here since the raw images are completely homogeneous with the 
other five. The full data analysis for NGC 7626 will be described 
in Harris et al. (2009, in preparation) along with four other giant 
ellipticals that have B— band images from ACS/WFC and /—band 
images from archival WFPC2 exposures. 



FWHM{GC) » FWHM {PSF). This is the 
case, for example, for HST imaging of clusters 
in M31 and other Local Group galaxies (e.g. 
iBarmbv et al.|[2007[) . 

• Partially resolved: The GC and stellar pro- 
files have comparable widths, FWHM{GC) ~ 
FWHM{PSF). This is the case for HST imag- 
i ng of many of the GCs in the Vir go galaxies 
(| Jordan et al.ll2005t iMieske et al."2006') at d = 16 
Mpc and for NGC 5128 (Gomez & Wo o dlevl [20071 : 
[McLaughlin et al.|[2008f ) at d = 4 Mpc, for exam- 
ple. 

• Marginally resolved: The GC is barely large 
enough to be measured relative to the PSF, 
FWHM{GC) - 0.1 - 0.3 FWHM{PSF). The 
GCs analyzed in the present paper fall into this 
category. 

• Unresolved: FWHM{GC) < 0.1 FWHM{PSF), 
for all practical purposes making it starlike and 
allowing normal PSF-fitting photometry. GCs 
i n much more di stant galaxies such as Coma 
(jHarris et al.ll2009( l. even with HST resolution, fall 
into this category. 

The mean half-light radius of classic GCs is typically 
about 3 pc, and 90% of them have < 6 parsecs (Fig- 
ure [1]). Thus at d ~ 40 Mpc, they will subtend a typ- 
ical angular diameter ~ 0.03", which equals about 30% 
of the 0.1" FWHM of stars on the HST ACS camera. 
This comparison means that we should expect the GCs in 
these target galaxies to fall in the marginally resolved or 
(for the smallest ones) unresolved regime. Although the 
differential aperture-size corrections to their photometry 
can be expected to be small, the purpose of the present 
paper is to determine these corrections and to construct 
a new database of photometry within which any MMR 
can be more securely determined. 

3. PHOTOMETRIC TECHNIQUES 
3.1. Defining the PSF 

Each target galaxy discussed in this paper is sur- 
rounded by some thousands of barely-nonstellar globu- 
lar clusters plus a variety of 'field' contaminants (a rel- 
atively few foreground stars plus numerous faint, small 
background galaxies). As will be seen below, the GC 
populations clearly outnumber all other types of objects. 
The first task for the photometry is to derive individual 
profile-size measurements of each detected object. PSF- 
convolvcd model GC profiles must be matched to each 
individual object and the optimum fit selected. Several 
codes have been developed and used in t he litera t ure to 
do this, including the ISHAPE code of [LarsenI (|1999D 
and iLarsen et al.l (|2001h which has been used for young 
and old star clusters in a wi de range of d is tant galax- 
ies; the KINGPHOT code of IJordan et aH HqoI) used 
pr incipally for the GCs in Virgo galaxies; and the code 
of iMcLaughhn et all (|2008l ) which can fit any one of five 
standard GC models and has been used for clusters in 
M31 and NGC 5128. 

The fitting code adopted for this work is ISHAPE. This 
code gives a fast, well tested, effective procedure for fit- 
ting marginally resolved GCs like these, where we are 
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primarily interested in determining the cluster size re// 
and cannot expect to obtain precise solutions for more 
subtle structural parameters such as ellipticity, orienta- 
tion, or core radius. These latter quantities become dra- 
matically easier to measure for star clusters in nearby 
galaxies, such as for Virgo, NGC 5128, or the Local 
Group members where the degree of resolution is much 
higher. ISHAPE and its parent package BAOLAB also 
contain highly effective tools for testing the results, in- 
cluding profile-subtracted images, and the ability to gen- 
erate simulated populations of PSF-convolved objects of 
arbitrary size and ellipticity. 

It is possible to construct PSFs for HST images 
through TinyTim modelling, but difhculties arise for 
appropriate choic e s of a diffusion ke rnel (see, e.g. 
iSpitler et all l2006t iGeorgiev etldl 120081 for recent ex- 
tensive discussions in similar studies) . To avoid any such 
uncertainties in the modelling chain I have employed 
purely empirical PSFs constructed for each individual 
frame. Fortunately, the true FWHM of the stellar PSF 
on the ACS/WFC frames is known beforehand to fall 
in the range 0.09" — 0.10" on series of images that have 
been properly coadded. This a priori constraint helps to 
select out candidate starlike objects. 

All exposures for each target field and in each filter 
were registered and comb ined through mult idrizzle, as 
described in more detail in [Harris et al.l (120061 ). To build 
a prelimi nary catalog of detected objects, the SExtrac- 
tor code ()Bertin fc Arnoutill996D was used with a gener- 
ously low threshold, and from this, a diagnostic graph of 
object half-light radius versus total magnitude was 
plotted. An example for the NGC 4696 field is shown in 
FigurelD The next stage was an iterative process of iden- 
tifying genuine starlike, isolated objects across the image 
from which the PSF could be defined. Objects selected as 
candidate PSF stars were those within the region marked 
out in the Figure, falling along the lower envelope of 
the ri/2 distribution. Moderately bright, unsaturated 
stars fall within this region, whereas saturated stars fall 
along the upturn region at the bright end of the dis- 
tribution. The candidates were then inspected through 
IRAF/imexamine, their radial profiles and FWHMs mea- 
sured, and ones with FWHM falling in the range 1.8 — 2.0 
px (that is, 0.09" — 0.10") were selected. This subsample 
was put into the IRAF/daophot/psf routine, inspected 
at finer detail, and any with remaining excess nonstellar 
wings or faint neighbors rejected. The final numbers of 
objects used in each field, along with the final FWHM of 
the adopted point-spread function, are listed in Table 2. 
In all cases, the stars used to define the final PSF were 
widely distributed in location, and thus the PSF repre- 
sents an average across the field (though no differences 
in the PSF shape or width at the level of ±0.1 px as a 
function of position were found; see Harris et al. 2006 
for additional discussion). 

Other characteristics of the raw image data including 
the internal photometric precision and complete n ess ar e 
described in much more detail in iHarris et al] (|2006f ). 
Briefly, the exposures for each galaxy were originally de- 
signed to reach an intrinsic limiting magnitude Mj ~ —8 
at 50% detection completeness, so that the turnover 
point (peak frequency) of the globular cluster luminos- 
ity function (GCLF) at Mj ~ —8.4 would be reached at 
almost 100% completeness. The data then thoroughly 



sample about 4 magnitudes of the GCLF, essentially its 
entire bright half. The internal photometric measure- 
ment uncertainty is smaller than a{B — I) = ±0.10 mag 
for Ml < —8.5 as well. This color spread is approxi- 
mately the intrinsic Gaussian width of the blue sequence 
and is smaller than the width of the red sequence, allow- 
ing the bimodal sequences to be well resolved over the 
full magnitude range of interest here. 

3.2. Measuring Effective Radius and Aperture 
Corrections 

As mentioned above, the single most important struc- 
tural parameter we are after is the effective or half-light 
radius of the GCs around these giant ellipticals. This 
quantity is what determines the amount of "aperture cor- 
rection" to apply to the measured magnitudes. In addi- 
tion, the cluster scale sizes are interesting in their own 
right as possible functions of galactocentric position and 
GC metallicity and will be discussed separately below. 

Within ISHAPE, r^ff is found by fitting the PSF- 
convolved cluster profile to each object in the image, then 
varying the assumed r^f / till a best fit is achieved. The 
code also returns a (S/N) ratio, which depends mainly on 
object brightness and can be used to assess the quality of 
the result. For any objects judged to be starlike, a value 
Te// = is returned. There are other structural param- 
eters that can be solved for in principle, including the 
ellipticity e, orientation angle 9, and central concentra- 
tion ratio, but these are easy to determine only for much 
more well resolved objects than we work with here. 

ISHAPE uses a lOx-subsampled version of the 
IRAF/daophot PSF as a basis for its image convolution. 
For the intrinsic GC profile, i t can use a variety of an- 
alytic functions; I adopt the iKingl (|1962| ) model since 
that model was defined originally to match real GCs 
and has been thoroughl y tested in other star-cluster pro- 
grams of this type (e.e. lLarsenlll999l; iL arseneLaLlllQQl 
ISpitler et al.l[200i IGeorgiev et al.l[2008h . The GCs in the 
target galaxies here are not nearly well enough resolved 
to allow ISHAPE to solve independently for the clus- 
ter concentration ratio c = log(rt/rc) where rt, Tc are the 
tidal and core radii, so a fiducial c = 1.5 ("KingSO" in the 
internal notation of ISHAPE) w as a dopted for all the so - 
lutions, foUowing lLarsenl (jl999f ) and lLarsen et al.l (|2001h . 
This c— value matches the mean King-model concentra- 
tion for the Milky Way, M31 and NGC 5128 globular 
clusters, among others ( Harrisl ll996l : iBarmbv etTalibOOTj : 
iMcLaughlin et all [2008) . In Fig. [11 the c-distribution 
for the Milky Way clusters is shown; 83% of the non- 
core-collapsed clusters are in the range 1 < c < 2. In 
addition, ISHAPE was run with the assumption e = 0, 
because the vast majority of known globular clusters have 
quite modest projected ellipticities, mo stly e < 0.1 (e.g. 
IHarris et al] 120021 : iBarmbv et al]|2007[) . Through a set 
of simulations to be described below, the sensitivity and 
accuracy of the r^ff solution to the object's brightness 
and assumed c, e were tested. In practice, ISHAPE 
returns the best-fit FWHM for the GC model profile, 
which must then be converted to an actual effective or 
half-light radius through a multiplicative factor that de- 
pends on c; for the specific case of the King30 model, 
r^ff =rh^ lASFWHM. 

In all cases except the few objects that happened to be 
affected by extreme crowding (neighbors within about 
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3 px), the residual maps plainly showed good fits to 
the PSF-convolved model. Some sample ISHAPE fits 
to three different objects are shown in Figure [31 

Knowing the value of r?i, a self-consistent procedure for 
defining the total magnitude of the object must be con- 
structed. The procedure adopted here was to start with 
a fixed-aperture magnitude for the object using rgp = 2. 5 
px, with the standard tools in daophot (jStetsonlll987f ). 
This fiducial 2.5-px radius is slightly larger than the stel- 
lar FWHM of ~ 2 pc but yet small enough to avoid 
crowding issues. This magnitude was then corrected to 
"large" aperture with a curve of growth (COG) of en- 
closed light versus aperture size; for the purposes of this 
paper, I adopt r^ax = 20 px = 1", at which the curves-of- 
growth for starlike and marginally resolved objects have 
accurately converged. The key point is that each value 
of effective radius has its own curve of growth, with more 
extended objects having relatively more of their light at 
larger radius. 

A complete grid of COGs generated through the 
ISHAPE simulation tools, for values of FWHM ranging 
from to 3 pixels, is shown in Figured This range cor- 
reponds to Te// ~ — 30 pc linear radius, a generous cov- 
erage of the sizes of known GCs and even UCDs. More 
realistically, we would expect the GCs in these galaxies 
to have radii ^ 6 pc = 0.6 px if they at all resemble 
known GCs. As can be seen from the Figure, this pri- 
mary range of interest is marked by the solid vertical 
bar in each panel. The figure shows that the COGs 
converge extremely well to the same total magnitude 
at the large (20-px) radius and the aperture corrections 
are small. It can readily be seen that this procedure 
is equivalent to making the 2.5-px aperture magnitude 
brighter by the amount needed to shift it onto the star- 
like curve{reff = 0); then, the correction from 2.5 px to 
20 px aperture radius is made along the starlike curve, 
which is the upper envelope of the grid of tracks in the 
Figure. 

Figure [5] shows the magnitude correction Am(2.5px — 
20px) as a function of object FWHM for the two filters. 
Again, the regime of interest here is almost entirely at 
the lower end (r^ff ^ 0.6px) where the aperture correc- 
tions are quite small. In practice, for about a third of the 
objects in the measured sample, the sizes turn out to be 
indistinguishable from stars, though most of these star- 
like objects are also very faint and of no interest in the 
following discussion. In addition, the correction curves 
in B and / are nearly parallel, so that the net aperture 
correction to the color {B — I) is very small. 

The images in the /—band go slightly deeper than in 
the B— band for objects with GC-like intermediate col- 
ors, so that the S/N of the ISHAPE solutions is higher in 
/. This comparison is illustrated in Figure[6]for the NGC 
4696 field, in which {S/N)b ^ OA{S/N)i. The results 
for all the other fields are quite similar. The COG appro- 
priate for each filter was used to define the total B and / 
magnitudes for each object. Finally, the S'/A^— weighted 
mean of reff{I) and reff{B) was adopted as the fidu- 
cial effective radius. The color index of each GC follows 
directly from {B - I) = Btot - hot- 

3.3. Tests of the GC Size Measurement 
3.3.1. Sensitivity to Central Concentration 



Although the King30 model accurately represents the 
mean for real GCs, the individual rt/rc ratios can differ 
between GCs by factors of 3 or more in either direction, 
or ±0.5 in c (see Fig.[T]). To test the sensitivity of the so- 
lutions to the adopted c, I repeated the solutions on sev- 
eral of the galaxy fields adopting much higher and lower 
c— values: in the ISHAPE notation, these were KinglO 
(c = 1.0), Kingl5 (c = 1.2), and KinglOO (c = 2.0). 
Sample differences for two of the fields are shown in Fig- 
ure [71 These tests show that for clusters with intrin- 
sically large radii r^ > 0.05" — 1 px, the choice of c 
is important and systematic differences in the deduced 
cluster size occur that grow with r^. However, for these 
rare big objects c can be allowed to be a free parame- 
ter within ISHAPE and explicitly solved for. At smaller 
radii {r^ < 0.05") the systematic differences are less than 
0.01". Said differently, for these tiny objects c cannot be 
solved for as a free parameter, but making the wrong as- 
sumption for c leads to only very small systematic errors 
in the deduced size relative to the correct value: if the 
assumed c is too small, r^ is very slightly overestimated, 
and if c is too large, r^, i s slightly underestimated (see 
also iGeorgiev et al.l 120081 . for similar tests and conclu- 
sions). The numerical reason why rh is rather insensitive 
to c in this small-object regime is simply that the con- 
version factor /(c) — (FWHM/rh) depends on c itself, 
and so an initial overestimate of the FWHM is mostly 
compensated for by an opposite change in /(c) to reach 
an equally good fit. In this "marginally resolved" regime, 
the object-to-object rms scatter and systematic shift to- 
gether add up to an uncertainty cr{rh) — ±0.003" due 
solely to the intrinsic uncertainty in c. 

3.3.2. Sensitivity to Cluster Ellipticity 

The assumption e = is clearly an idealized case. If 
a cluster is actually elliptical, how much of an effect will 
that assumption have on the deduced r/i? Figure [51 shows 
the results of an extensive set of simulations. Here, arti- 
ficial clusters of intrinsic shape e = {1 — b/a) = 0.2 over 
a wide range of brightness and a wide range of sizes were 
convolved with the PSF, put onto an image with realis- 
tic background noise, then remeasured with ISHAPE. It 
should be noted that the adopted value e = 0.2, though 
modest, is already much larger than the ellipticities of 
most real GCs, so that any bias that shows up in these 
tests will be an upper limit to the average for real clus- 
ters. 

In Figure [51 the measured FWHM from ISHAPE is 
plotted as a function of S/N (essentially the brightness 
of the object). In the left panel, ISHAPE is allowed 
to solve for e = 1 — {b/a) as a free parameter, and in 
the right panel, the code is required to assume e = 0. 
The four different rows show input objects from big to 
small, with semimajor axis from 1.0 down to 0.1 px (that 
is, from 50% down to 5% of the PSF width). In both 
cases, the true value of the FWHM is accurately returned 
for S/N > 50 and for intrinsic sizes bigger than about 
10% of the PSF FWHM. For fainter and smaller objects 
than these limits, the solution for the object size becomes 
indeterminate. 

Additional tests of the accuracy of the solutions are 
shown in Figure[51 Here, the results from the same simu- 
lations are shown for the measured axial ratio b/a and for 
the orientation angle 6. For these simulations the "right" 
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answers (the input values) are {b/a — 0.8, 9 — 45°). The 
results of these tests show once again that for the brighter 
objects with S/N > 50 and intrinsic sizes FWHM > 0.2 
px (about 10% of the PSF size), accurate answers are 
returned. For the faintest objects, however, the spread 
of solutions for FWHM, (b/a), and especially 9 rapidly 
blows up. For the smallest objects as well (FWHM less 
than 10% of the PSF width), the solutions for object 
shape and orientation are unreliable. 

The random measurement uncertainty for = 
1.48 FWHM in the regime S/N > 50 and rh < 0.5 
px is, from the simulations, a{rh) = 0.0025". No sys- 
tematic errors in the solutions for the object size, shape, 
and orientation occur until FWHM < 0.1 px, at which 
stage the object can be treated as virtually starlike. 

3.3.3. Intrinsic Scatter and Total Uncertainty in Size 
Measurements 

A convenient summary of the range of object sizes that 
ISHAPE can be expected to solve for is shown in Fig- 
ure [ini This is the result of a separate set of simula- 
tions where the input objects all have {S/N) > 50 and 
a range of sizes extending from "starlike" (FWHM=0) 
up to FWHM=5 px. Just as in the previous tests, the 
model cluster has a King30 profile and it is convolved 
with a PSF with width 2.0 px. Each point shows the av- 
erage of many dozens of simulated objects, with the error 
bar representing the standard deviation of the individual 
scatters around the means. The graph shows that for 
these relatively bright objects, the measured sizes are 
systematically accurate down t o 10% of the P SF size, in 
agreement with the findings of iLarseiil ()1999( ). 

The consistent message from all these tests is that we 
can accurately recover the true size and shape of the ob- 
ject if it is reasonably bright and large enough to satisfy 
FWHM > 0.1 FWHM(PSF). For objects in this "well 
conditioned" range, the total uncertainty in our size mea- 
surements can be summarized as follows. 

(a) The random uncertainty in the object size due only 
to the ISHAPE fitting procedure, as noted above and 
shown from the scatter in Figs. [8l and [TOl translates to a 
mean of ±0.0025". 

(b) The additional uncertainty due to the possible range 
of c— values (Fig.[7|) averages ±0.003". 

(c) Finally, there is an external uncertainty due to the 
PSF size itself. The numbers listed in Table [5] show that 
for each galaxy field and filter, the actual PSF width 
is internally uncertain by typically ±0.0033" in B and 
±0.0018" in /. These translate to an ISHAPE fitting 
uncertainty in of about ±0.005". 

Summing these three independent effects in quadra- 
ture, we obtain a net random uncertainty per object of 
a{rh) = ±0.006" or ±0.12 pixel. For a galaxy at d ~ 40 
Mpc, this converts to a true uncertainty in scale size 
of typically a{rh) — ±1.1 parsecs due solely to the GC 
profile-fitting procedure. 

A comparison of the measured object sizes, done inde- 
pendently on the B and / images, is shown in Figure ITT] 
for each of the six galaxy fields. In most cases, the mu- 
tual agreement between the two filters is good and well 
within the expected differences due to PSF-size uncer- 
tainty. The one exception is for NGC 1407, where rfi{B) 
is systematically bigger than rh{I) by ~ 0.2 px, an offset 
almost twice as large as the expected uncertainty due to 



the PSF sizes mentioned above. The first source of this 
anomaly to suspect would be that the PSF size in this 
particular field is too big in /, or too small in B, or some 
combination of these; however, the PSFs in both filters 
have very much the same intrinsic widths as in the other 
five fields (Table 2) and were derived with exactly the 
same careful iterative procedure. Further tests have not 
revealed any clear solution and I have chosen to leave 
the results as they are, at very least to ensure that the 
final procedures for all the fields are homogeneous. The 
net effect on the mean cluster sizes for NGC 1407 is not 
as large as this comparison might indicate, because the 
final rfi value is the weighted average from the two filters 
and the / filter carries most of the weight (see Fig. 

Fig. [11] also allows us to make one last test of the inter- 
nal measurement uncertainties. The FWHM measure- 
ments are done separately for the same objects on the 
B and / frames, with independently determined point- 
spread functions, so the scatter in the B, I size measure- 
ments provides a reasonable estimate of the measurement 
uncertainty. The rms scatter is closely similar over all 
six fields, and averages a{rh) = 0.16 px, equivalent to 
±0.008". Accounting for the fact that the final rh values 
are the average of B and / then reduces the uncertainty 
to ±0.006", in close agreement with the estimate given 
above. 

4. THE COLOR-MAGNITUDE DISTRIBUTION 

With the object size measurement and aperture pho- 
tometry complete, we now have a database of (_B, /) pho- 
tometry where the magnitudes of each individual object 
were determined starting from small-aperture photom- 
etry and then fully size-corrected to total magnitudes. 
The resulting distributions in the raw color-magnitude 
plane (/, B ~ I) are shown in Figure [T^ 

In each galaxy the bimodal GC sequences are evident, 
but the first question to deal with is to ask how these new, 
recalibrated data match up with the previous ph otome- 
try based primarily on PSF-fitting (jHarris eFaDIlooe). 
A graphical comparison is shown in Figure fT3l The huge 
statistical weight of the combined data, compared with 
any other GCS study, can be realized by noting that this 
diagram contains 15000 objects that were measured by 
both psf and corrected-aperture methods, most of them 
globular clusters. 

Partially anticipating the results of the later discus- 
sion, we use the mean lines for the blue and red GC se- 
quences in Fig. [13] (calculated as described later in Sec- 
tion 16. ip to note that both the old and new datasets 
define closely similar bimodal sequences. In other words, 
completely reworking the data according to a technically 
improved procedure has yielded no first-order change. 
The major reason for this similarity is that the cluster 
sizes are small enough that the individual aperture-size 
corrections relative to the fiducial starlike curve (Fig. [4]) 
are modest. When they are differenced between B and 
/, the aperture-size corrections to the color indices are 
even less. 

To a finer level of detail, comparison of the two datasets 
shows that the final colors in the present size-corrected 
data are a few hundredths of a magnitude bluer than 
in the older data, particularly in the bright range Mi < 
— 10. As will be seen in the next sections, this offset is not 
due to some particularly strong size/luminosity relation. 
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because the characteristic sizes of the GCs have large 
scatter at any luminosity, and they are all individually 
accounted for. The net result is that the deduced MMR 
on the blue sequence (the slope of the correlation between 
luminosity and color) is still present but at slightly lower 
amplitude than in the earlier data. These outcomes will 
be discussed in more detail below. First, however, we 
look more thoroughly at the distribution of sizes of the 
clusters themselves, which have several points of interest 
in their own right and are an important preamble to the 
features of the MMR. 

5. ANALYSIS OF THE CLUSTER SIZES 

The characteristic (half-light) radius of a GC is 
an astrophysically valuable and frequently used quan- 
tity because it is expected to stay nearly invari- 
ant with internal dynamical evolutio n well after the 
initia l cluster formation period (e.g. ISpitzer fc ThuanI 
19721; lAarseth & Heggid Il998t iBaumgardt et al.l 120021 : 
Trcnti et al. 2007,). and so the measured rh should 
represent an intrinsic structural property built in at 
an early time. It has long been realized that 
rfi shows large cluster-to-cluster differences but on 
the average it increases systematically with (three- 
dime nsional) Galactocentric dis tance Vgc in the Milky 
Way ([van den Bergh et al.lll99l[ ) following a rough scal- 
ing rule rfi ~ ''^g'c- This trend is generally interpreted 
as meaning that proto-GCs in the process of formation 
will take up systematically larger scale sizes if form- 
ing within shallower surrounding potentials. Similar, 
though weak, correlations in the observed over a large 
range in projected (two-dimensional) galactocentric dis- 
tance Rgc have been found in oth er galaxies such as the 
edge-on disk galaxies NGC 5866 ('C antiello et~al] l2007| ) 
and NGC 4594 (Soitlcr ct al. 2006), and the nearby gE 
NGC 5128 (iGomez fc Woodlevll2007D . typically scahng as 
rh Rg'^ . A similar trend can also be found within 
the c omposite data for the Virgo galaxies (jJordan et al.l 
|2005() , though the data in this case cover a much smaller 
range in Rgc- 

In addition, differences in the mean r^ be- 
tween red (metal-rich) clusters and blue (metal- 
poor) clusters at the sa me galactoce n tric di stance 
have been noted (see iKund u et al.' '1998', '19991; 
Larsen et al.ll200lH Jordan et al.ll 2GG5: Spitlcr ct al. 200S; 
Gomez fc WoodlevI l2007l for a samphng of the key re- 
sults). The effect claimed in these papers is for the 
metal-richer clusters to be about 20% smaller than the 
blue ones. Although considerable scatter in r^ is found in 
each subgroup, this mean difference has persisted as the 
size and quality of the available measurements have in- 
creased. This size difference has been s uggested to be the 
resul t of a geometric projection effect ()Larsen fc Brodid 
l2003f ) whereby the metal-richer clusters lie in a more 
centrally concentrated distribution within the galaxy, so 
that a higher proportion of them are at smaller Rgc and 
thus physica Uy smal l er bec ause of the rh/Rgc correlation. 
Alternately, iJordanI ()2004D suggests that the difference 
is an intrinsic function of metallicity and the result of 
metallicity-dependent stellar evolution times. 

Data that cover a large range in Rgc and that are 
based on large statistical samples of clusters are needed 
to sort out the interpretations. If the difference is 



due to a project i on effe ct along the lines suggested by 
iLarsen fc Brodid ([2003), then it should become much 
smaller at larger distances outward into the halo. 

5.1. The GC Size Distribution 

The database of rt values from the six giant ellipticals 
in this paper clearly does not have highly reliable indi- 
vidual measurements of cluster size. But it has the com- 
pensating advantages of large sample size and coverage 
of a large range in galactocentric distance (Rgc/Reff)- 
The measurements for our six fields are displayed indi- 
vidually in Figure [TH Here, the half-light sizes are shown 
for all objects with r/j > 1 pc regardless of their internal 
quality (that is, not selected by (S/N); the particular 
cutoff of 1 pc will be justified below). In every graph, 
the bimodal nature of the color distribution is evident, 
as well as a consistent pattern for the blue GCs to extend 
up to larger sizes than do the red GCs.^ 

If the GCs in these giant ellipticals are basically the 
same kind of star cluster as the classic Milky Way clus- 
ters, then we would expect their size distribution as a 
whole to resemble the Milky Way's. In Figure [151 the 
combined results for all six fields are shown in histogram 
form. Here, the database is sampled in a different way 
than in the preceding Figure: objects of all colors are in- 
cluded regardless of their r^ value, but only the highest- 
quality measurements {S/N > 50) are included. This 
selection leaves 3330 objects over all half-light radii. For 
Tfi ^ 2 pc (the shaded region, equivalent to 0.01" at a 
distance of 40 Mpc), we approach the limit of measure- 
ment of the profile fitting code and any smaller values for 
the intrinsic size returned by ISHAPE become very un- 
certain. The typical internal uncertainty per object, from 
the arguments in the previous section, is ±1 pc, shown 
as the error bar at upper right. In other words, for the 
extreme case r/j < 1 pc, GCs cannot be distinguished 
from stars even under the most optimistic reading of the 
data. However, it should be noted that GCs this com- 
pact physically appear to be quite rare (there are only 
two this small in the Milky Way). 

While keeping in mind the measurement limits men- 
tioned above, the total histogram shown in Fig. [15] clearly 
has two major components: the dominant GC population 
peaking at r^ — 2 pc; then a dip near r^ = 1 pc and a 
secondary peak at r;i = which should consist of all the 
starlike or near-starlike field contaminants plus the most 
compact GCs. As a further test, ISHAPE was run on 
a simulated population of 5300 starlike objects {rh = 0) 
over a wide range of brightnesses mimicking the range of 
the real data. The size measurements of these stars are 
shown in Fig. [15] as the green histogram. The width of 
this histogram should give a reasonable estimate of the 
internal precision of the profile- fitting results. The vast 
majority of these simulated stars fall within rh < 1 pc 
and confirm the empirical internal error estimate of ±1 
pc deduced above. 

These various tests of the data show that the total 
spread of the major GC component that continues on 
upward past rh 2 pc is too large to be explained simply 

^ The mean values for NGC 1407 are noticeably smaller, a 
likely result of the discrepancy between the PSF sizes in B and I 
as discussed above. However, this galaxy provides less than 10% 
of the total GC population in the study, and whether or not it is 
excluded has no important effect on the subsequent discussion. 
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by measurement scatter. It is completely consistent with 
the interpretation that we have genuinely resolved the 
intrinsic size distribution of the globular clusters around 
these galaxies. A final key test is to compare it with the 
Milky Way GC sizes. In Fig.[T5]the Milky Way data from 
[Harris (,1996) are shown in the dotted histogram for the 
66 clusters with luminosities My < —7.3 (i.e., brighter 
than the turnover point of the GC luminosity function), 
and no further from the Galactic center than 20 kpc. 
Selecting the Milky Way sample this way makes it closely 
comparable to the GC sample that we have in the giant 
elliptical data. The overall shapes of the two datasets 
resemble each other remarkably well, particularly their 
extension to large radii. 

A second-order difference is that the Milky Way Vh 
distribution peaks at a point about 0.5 pc bigger than the 
peak for the giant elliptical data. This offset corresponds 
to an angular size difference of Ar?i ~ 0.0025" at d ~ 40 
Mpc and is entirely within the internal uncertainties of 
the measurements in this paper. // the offset is physically 
real, however, it is likely to be because the gE's have 
a much higher proportion of red (metal-rich) GCs than 
does the Milky Way, and as will be seen below, the red 
GCs are systematically smaller than the blue GCs by 
~ 0.5 pc (see Table [3] for the median values in the Milky 
Way). 

5.2. Correlations of GC Size with Other Factors 

Having presented the overall size distribution, we now 
investigate correlations with the three properties that 
have been suggested in the previous literature to influ- 
ence GC scale size: luminosity, metallicity, and galacto- 
centric distance. Differences in GC scale size with metal- 
licity have drawn the most attention in recent work, 
primarily from t he HST-based photome try of GCSs in 
nearby galaxies (iKundu et al.l IT999: Lar sen et al.l I200l1 : 
I Jordan et al.l 120051 among several others). These same 
studies have not, however, been used to track size cor- 
relations with either luminosity or spatial location very 
effectively, because (a) as will be seen below, the size 
difference with cluster luminosity becomes obvious only 
at the top end of the luminosity range where GCs are 
rare; and (b) these previous studies were almost entirely 
restricted to the inner-halo clusters within about 1.5 ef- 
fective radii of the galaxy spheroid light, so that only a 
very limited run of spatial location could be traced. The 
database in this study has the advantage that all three 
factors can be traced at least to some degree. 

In Figure 1161 the r^ distributions are shown now for 
the 2454 GCs bigger than 1 pc and with (S/N) > 50. 
The right panel shows the histograms separately for the 
metal-poor GCs (those with {B — /)o < 1.8) and the 
metal- rich ones {{B — I)o > 1.8). The expectations from 
the previous literature are borne out here, that the lower- 
metallicity GCs are systematically bigger on average by 
~ 0.5 — 1 pc, or roughly 15 — 20%. It is highly unlikely 
that this difference is just an artifact of measurement, 
because all the objects in the database, red and blue 
alike, were measured through both B and / filters and 
appropriately S'/A^— weighted. Both types of clusters also 
appear at all galactocentric distances and over the full 
range of GC luminosities allowed by the magnitude limits 
of the photometry. 

One caveat in these kinds of comparisons is that the 



sample has been deliberately truncated for < 1 pc. 
This should, however, have only a small effect on such 
diagnostic quantities as the median size of a subsample, 
because at least in the Milky Way virtually no clusters 
have such ultra-compact structures (Fig. [T5|) . 

The distribution of sizes versus cluster luminosity, 
shown separately for the blue and red GCs, is in FigurelTTl 
and the median sizes in half-magnitude bins are listed in 
Table |3l As in the previous figure, the 2454 objects with 
(S/N) > 50 and r^ > 1 pc are shown. For luminosities 
higher than Mj ~ — 11 (corresponding to about lO^L©) 
the clusters show a weak trend to become larger with 
increasing luminosity. For the lower-luminosity range 
Ml > —11, any trend of median cluster size with lu- 
minosity is modest at best and perhaps not significant. 
Both of these results agree wit h what has been fou nd in 
the Local Group galaxies fe.g. iBarmbv et ai]|2007i ). 

On average, the median GC size of the blue clusters 
is 15-20% bigger than the red clusters, at any luminos- 
ity. A standard Kolmogorov-Smirnov two-sample test 
also shows that their size distributions are different at 
> 99% significance. 

In the previous literature, correlations of size with 
galactocentric distance have also been suggested. Even 
within the Milky Way, which is the only galaxy for which 
we can explore that correlation in three dimensions, the 
trend is relatively modest, and when it is projected onto 
two dimensions for any other galaxy it can be expected 
to be even shallower. Nevertheless, the data in this pa- 
per can be used to gain some idea of any such correlation 
because the range of Rgc is quite large for each of the six 
galaxies. To normalize all the data to a common system 
and allow treatment of all the galaxies together, we ex- 
press the projected Rgc value for each cluster in units 
of Reff, the effective radius of the galaxy's spheroid as 
determined from its surface brightness profile (see Table 
[Hand the NED). 

A schematic representation of the combined effects of 
spatial location and metallicity together is shown in Fig- 
ure 1181 where size versus color is plotted in three dif- 
ferent spatial zones: the inner halo {Rgc < Rcff), the 
mid-halo {Reff < Rgc < '^Reff), and the outer halo 
{Rgc > '^Reff, extending out as far as 5Rcff). A more 
quantitative comparison is shown in Figure [T9| and Ta- 
ble [31 Here, the median values in bins of AR — 0.5 
are plotted separately for the two metallicity subgroups. 
These show a weak trend for both types of clusters to 
have systematically larger scale sizes farther from their 
galactic center. Furthermore, the slopes of both relations 
are quite similar: least-squares solutions for the medians 
give 

rh.med{pc) = (2.236±0.079)-f (0.203±0.030)(i?gc/i?e//) 

(1) 

for the blue clusters, and 

rh.med{pc) = (1.881±0.043) + (0.181±0.016)(i?gc/i?e//) 

(2) 

for the red clusters. Consistent with what was found 
earlier, the blue GCs are larger by 17% than the red 
GCs, and this difference does not change significantly 
with Rgc bin. 

A power-law for mu lation of the sarn e trend 
(Ivan den Bergh et al.l 119911: iSpitler et all l2006t 
iGomez fc WoodlevI l2007i ) has perhaps a stronger 
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physical motivation and yields an equally good repre- 
sentation of the data. The best-fit equations in this 
form (for the blue and red clusters respectively) are also 
shown in Fig. [19] and are given by 

; (pc) = 2.53{Rgc/ReffT '' , (3) 

(4) 



(pc) = 2.l5{Rgc/Reff) 



There are only two other galaxies beyond the Lo- 
cal Group for which GC-size measurements cover a 
similarly large run in s patial location: NGC 5128 
([Gomez fc Woodlevl 120071 ) , with data extending out to 
~ SReff, and NGC 4594 ([Spitler et a l. 2006), with data 
out to ^ 6i?e//. Both of these galaxies are relatively 
nearby and thus their GCs are better resolved, but the 
sample sizes are about one order of magnitude smaller 
than the one used here. Nevertheless, these two studies 
find power-law dependencies for rh{Rgc) with very much 
the same power-law slopes ~ 0.1 — 0.2 as found here. In 
short, all the current data suggest a consistent pattern 
for a shallow but steady outward increase in cluster scale 
size regardless of metallicity. 

5.3. Interpretation 

A brief summary of the findings for the cluster scale 
sizes in these six giant ellipticals can be made as follows: 

• The characteristic half-light radius depends weakly 
and nonlinearly on total luminosity: for the lower- 
luminosity range that contains the vast majority 
of the clusters, i.e. the range M/ > — 11 or i < 
1O^L0, the median for clusters of all metallicities 
remains largely independent of luminosity, whereas 
for the bright end M/ ^ — 11 the median GC size 
begins to increase steadily. For the extremely rare 
and very most luminous GCs at Mj ~ —12.5, the 
scale sizes are almost twice as large as the median 
for the lower- luminosity range (Fig. I17p . 

• At any luminosity, and at any galactocentric dis- 
tance, the metal-poor (blue) clusters are system- 
atically larger than the metal-rich (red) ones by 
(17 ±2) percent (Figs. [13 [H and[Tl|). 



R%}j (Figs. [18] and 



• Lastly, there is a consistent correlation of GC scale 
size with projected galactocentric distance, at any 
metallicity, equal to r/j med 
[TU). 

Given the limited material that was available from pre- 
vious work, the fact that all three correlations are compa- 
rably weak explains why it has been difficult to separate 
them from one another, or indeed to decide whether or 
not each one is even present. 

The systematic increase in scale radius with luminosity 
for the highest-mass clusters already has clear evidence 
in its favor from the Local Group galaxies fBarmb v et al.l 
l2007f ) and the brightest GCs and Ultra-Compact Dwarfs 
(UCDs) in the Virgo and Fornax ga l axies (iHasegan et al.l 
[20f)5tlEvstigneeva et al.l[2nf)7ll2008h . iSpitler et al.l (120061 ) 
also found a significant increase in mean size with lu- 
minosity for the GCs in the Sombrero galaxy, with a 
noticeable "onset" point near My — —9 above which 
the effective radii begin to increase most strongly. These 
studies show repeatedly that these compact, massive stel- 
lar systems have (r^) const for luminosities less than 



Ml — 11 = IO^Lq, and then gradually swing upward 
to a scaling (r/i) ~ Ly*" which bridges them to other ob- 
jects such as d E nuclei and eventually dw arf galaxies (see, 
e.g., Fig. 2 of lEvstigneeva et al] (|2008f ) for a useful and 
up-to-date composite diagram). Nevertheless, the cur- 
rent data show few luminous GCs that match the larger 
sizes associated with UCDs. The mean size/luni i nosit: 
relation for UCDs, taken from lEvstigneeva et al.l (|200' 
and converted to M/ with {V — I) = 1.0, is plotted as 
the dashed lines in Fig. [T7] Most of the GC data from 
these giant ellipticals lie well below the UCD curve even 
at the highest luminosities. In other words, most of the 
"classice" luminous GCs have compact structures that 
are smaller than UCDs. These GCs with Mj < -11.5 in- 
stead more closely matc h the lo wer envelope of the data 
discussed by lEvstign eeva et al.l ((2008) . This lower en- 
velope is populated mostly by a handful of dE,N nuclei. 
This comparison, perhaps, provides additional hints that 
the most luminous normal GCs may have had their ori- 
gins as the nuclei of dwarfs that were later stripped as 
they were accreted by the giant elliptical. 

It should be kept in mind, however, that in a normal 
GC luminosity distribution with turnover at M/ ~ —8.3 
and dispersion ~ 1.4 mag, only ~ 1% of the entire 
GC population lies above Mj < —11. For the lower- 
luminosity GCs that make up the vast majority, the as- 
sumption r/i ^ const continues to be a useful one. 

The present study was not especially tuned to mea- 
suring the sizes of larger UCD-like objects, so a separate 
analysis of them in these fields will have to be conducted. 
Work of this kind has begun and will appear in a later 
paper. 

In the model of lLarsen fc Brodid (|2003l ) . the mean size 
difference between metal-poor and metal-rich clusters is 
a geometric projection effect: if the intrinsic size distri- 
bution for all clusters is actually the same, and if they 
all follow the same r^irgc) trend, then the metal- rich 
ones will be smaller on average because their spatial dis- 
tribution is more centrally con centrated to the gal actic 
center. The later discussions of ' Spitler et al.l (l2006l ) and 
iGomez &: W oodlcv (2007) supported this interpretation. 
By contrast, Jordanj (j2004l ) presents an alternative model 
whereby the metallicity-based size difference is the end 
result of many Gy of dynamical evolution and mass seg- 
regation within the GC combined with stellar-evolution 
timescales that depend on metallicity. Yet another pos- 
sibility is that metallicity may have a more direct and 
much earlier effect on cluster size, in the sense that higher 
heavy-element abundance in the protocluster gas may en- 
courage more rapid cooling and cloud contraction before 
the stars form, allowing a metal-rich cluster to have a 
smaller scale size from the start. Detailed models along 
t hat line have yet to be p ursued. 

iLarsen fc Brodid ()2003f ) suggest that if the projection- 
effect model is correct, then the sizes of high- and low- 
metallicity clusters should become the same at large 
galactocentric distance well outside the galaxy core. Un- 
der these assumptions, their quantitative models show 
that the red-vs.-blue size difference should be greatest 
within one or two Reff- However, the data here (see 
especially Fig. [T9]) indicate that this is not the case: 
metal-rich GCs are systematically smaller than metal- 
poor ones by much the same factor at all locations, fa- 
voring the view that the metallicity-based difference is 
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intrinsic. At the same time, both types of clusters fol- 
low a size/distance dependence rn ^ Rgc with a similar 
power-law scaling exponent a, but just with different ze- 
ropoints. This effect strongly suggests that metallicity is 
not the whole story. The scale sizes of the clusters are 
at least partly determined by the depth of the tidal field 
that they are formed in, ending up with more extended 
structures in environments of shallower external poten- 
tial wells, such as in dwarf galaxies or the outer halos 
of large galaxies. Notably, in a recent survey of GCs in 
many nearby dwarf galaxies, Georgiev et al (2008, 2009) 
find a median = (3.2 ± 0.5) pc, and r^j ^ 5 pc in the 
LMC. Sizes this large are reached in our gE sample only 
for the outer regions at several R^f / . 

A final overview of this part of the discussion is that the 
scale sizes of globular clusters appear to be determined 
jointly and universally by their metallicity, spatial loca- 
tions, and (at very high luminosity) their mass: all three 
of these effects are physically real and need to be included 
in their conditions of formation. Deeper modelling of the 
combined effects of formation, space distribution, and 
dynamical evolution is far beyond the scope of this pa- 
per, but the sum total of all available data now appears 
to provide enough evidence to motivate such models. 

6. THE MASS/METALLICITY RELATION 

The final magnitudes and colors of the globular cluster 
populations in our six gE galaxies, as shown in Figs. 1121 
and [131 &re based on fully size-corrected magnitudes and 
can now be used to investigate the blue and red GC se- 
quences in more detail. The important questions, beyond 
the first-order one of establishing the basic bimodality of 
the GC color distributions, are to ask whether either se- 
quence shows a detectable mass/metallicity relation; and 
whether the slope or form of that relation might itself de- 
pend on cluster luminosity. Beyond that, we would like 
to test its behavior versus other external characteristics 
such as cluster size or galactocentric distance. 

6.1. Bimodal Fitting to the GC Sequences 

The approach to defining the MMR adopted here is to 
measure the mean color of each sequence i ndependently 
at many points along the luminosity range ([Harris et al.l 
[2006; Mieskc ct al. 2006). In some other previous studies 
JStradcr et al. 2006; Soitler et al.' '2006'; 'Cantiello^eLaD 
[2007; DcGraaff et al. 2007; Spitler ct al. 2008), mean Hn- 
ear relations in color versus magnitude were fit to each 
sequence and the best-fit slopes were quoted. If the GC 
sequence has the form A(color)/A(magnitude) = const, 
then this translates readily to a power-law form in heavy- 
element abundance versus GC luminosity, Z ^ L^, as 
long as the color index varies linearly with [Fe/H]. How- 
ever, the linear-fit approach imposes an uncomfortably 
strong assumption on the data. The raw data pro- 
vide no a priori guarantee that the mean GC colors 
will vary linearly with magnitude, and a forced linear 
fit will miss any nonlinear trend that might be present. 
Whether or not the MMR disappears below a certain 
mass range, for example, is important for key fea tures of 
the self-enrichment model ([Bailin &: Harrisl[2009| ). Fewer 
assumptions are forced on the material by finding a bi- 
modal fit to the distribution of GC colors in independent 
magnitude bins, and then inspecting the resulting mean 
points for any trends, linear or otherwise. 



The tradeoff is that if the MMR actually has a linear 
or near-linear form over the entire luminosity range, then 
its slope can be determined by the linear-fit method with 
a smaller number of objects than is required for the bi- 
modal fitting in independent bins. The approach I adopt 
here is that this linear-full-range assumption first needs 
to be tested independently on large, statistically unam- 
biguous datasets before being carried over to smaller 
galaxies. 

Another factor capable of biasing the deduced MMR is 
that of sample field contamination. In the present study, 
contamination by (mostly) faint background galaxies and 
a few foreground stars over a wide range of colors is 
nearly negligible over the brighter part of the GC pop- 
ulation but grows to significant levels for / > 24.5 (see 
Fig. [T2|) . The contamination on the blue side is typically 
larger and thus may leave the artificial impression of a 
too-steep MMR slope if a linear fit is imposed. A no- 
ta ble exception from previous literature is the database 
of lMieske et al.l ([2006( 1 where the GCs in the Virgo galax- 
ies were individually selected and the sample has low con- 
tamination. In general, however, bimodal fitting to the 
color distribution will be less affected by contamination 
than the linear-fit method since it is designed to pick out 
the peak color of each component regardless of the scat- 
ter of outlying points. The level of field contamination 
for the present data is very small for our main luminosity 
range of interest (M/ < —8.5, brighter than the GCLF 
turnover) but the issue will be discussed in greater detail 
in Section [6721 below. 

Th e RMIX multicomponent fitting code ([MacDonaldl 
l2007f) is used here to determine the bimodal fits. The 
code and its appl ication for thi s prob lem are described 
in more detail in IWehner et al.l ([20081 ). In short, start- 
ing from either the individual datapoints or a previously 
defined histogram, RMIX finds the best-fit multimodal 
description of the color distribution. The user-specified 
input consists of the number of modes to be solved for; 
initial guesses at their means and standard deviations; 
and the basic fitting function to be used. The usual 
Gaussian model is most common, but half a dozen other 
well known functional forms are available within RMIX. 
The mean, standard deviation, and proportion of each 
component can either be fixed by the user or solved for 
in any combination. 

The solutions for the six individual galaxies are sum- 
marized in Table [5] and shown in Figure [20l Here, mag- 
nitude bins Al of at least 0.5 mag were used, but with 
the proviso that a minimum of roughly 100 objects fell 
into each bin. To keep the fewest possible external con- 
straints on the solutions, RMIX was required only to 
assume a bimodal Gaussian distribution in (B — /)o and 
then was left free to solve for the mean {B — I) col- 
ors fj,i{blue), fj,2{red), the standard deviations of each 
mode tTi , (72 , and their proportions of the total popu- 
lation pi + P2 = 1. The input data included all objects 
within the color range 1.5 < [B — I) < 2.5, an interval 
chosen to reject the largest part of the outlying field con- 
taminants while safely including both GC components. 
In a very few cases (individual bins here and there) it was 
necessary to constrain the fit in order to get convergence, 
usually by fixing the proportions of one or the other of 
the ct's; these few cases are marked in Table[5]with inter- 
nal uncertainties of zero. Overall, however, this solution 
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technique differs very noticeably from most previous ones 
in the literature because the large GC populations allow 
these minimal constraints on the fits. 

The maximum statistical weight can be obtained by 
combining all six systems into a single dataset, where 
everything is translated to the common scale of absolute 
magnitude and intrinsic color M/, (B — /)o. Biniodal 
RMIX fits for the combined data are listed in Table [H 
These are also plotted in Figure [5T1 as are the mean 
points for the six individual systems. 

Assuming that the photometry itself is systematically 
correct, the clear conclusions from these bimodal fitting 
results are that 

(a) A nonzero relation between luminosity and color ex- 
ists along the blue-GC sequence particularly; 

(b) The trend is nonlinear in form; and 

(c) Along the red sequence, any such relation is gentler 
or absent. 

I show in Figure [51] the CMDs for those objects with 
positive, nonzero size measurements from the ISHAPE 
fits (that is, objects with r^, > 0), and then separately 
for all the rest of the objects that were either found to be 
starlike (r^, = 0) or too faint for successful ISHAPE fit- 
ting. The vast majority of objects in the second category 
populate the area of the CMD fainter than Mj = —8.5 
and thus do not affect anything but the faintest bin of 
the RMIX bimodal sequence fits. Removing the starlike 
objects entirely (left panel of Fig. [H]) does not change 
the best-fit red and blue GC sequences. 

6.2. Background Field Contamination 

Although the CCS populations clearly dominate the 
color-magnitude diagrams in the six galaxies studied 
here, it is useful to gauge the actual level of field con- 
tamination. The raw data in this HST/ACS program 
contain no true "control field" exposures pointed at adja- 
cent blank fields. However, three of the galaxies imaged 
in the original program turned out to have small GC 
populations, and the outskirts of their fields can be used 
for an estimate of the field populations. T hese are NGC 
5322, 5557, and 7049 (iHarris et al.][200l . The advan- 
tages of using these fields are that they were taken with 
the same filter combinations and to the same absolute- 
magnitude limits, and the galaxies themselves have sim- 
ilar distances, foreground reddenings, and Galactic lati- 
tudes to the other six. 

The fields of these other three target galaxies were all 
remeasured with the same procedures through daophot 
and ISHAPE as described above. Inspection of the color- 
magnitude diagrams showed that GC populations around 
them arc definitely present in the inner regions, but die 
off with increasing radius until for Rgc > 125" the GC 
numbers begin to fall below the general density of field 
objects that are widely distributed in color. Thus to 
generate a composite "control field" I have combined the 
regions beyond Rgc — 125", which add up to a total area 
11.2 arcmin^, conveniently equal to the area of a single 
ACS/WFC field. 

The CMD for this synthetic control field is shown in 
Figure [531 Here the raw (/, B — I) data have been con- 
verted to (M/, (B — /)o) assuming (m — M)j ~ 33.1 
and Eb~i — 0.05, which are the averages for NGC 5322, 
5557, and 7049. The data for "starlike" objects {ru = 0) 
scatter uniformly across the CMD and should represent 



genuine field objects. By contrast, the data for the ob- 
jects with rh > should include some genuine GCs that 
belong to these three galaxies (at a distance of 40 Mpc, 
Rgc = 125" corresponds to a linear projected distance of 
only 25 kpc, which is still well within the bounds of typ- 
ical gE halos), along with a variety of small background 
galaxies. Traces of the red and blue GC sequences can 
indeed be seen in the middle color range. In other words, 
the composite field data can represent only an approxi- 
mate upper limit to the true level of field contamination 
in the six main target galaxies. 

The distributions in color for the background data are 
shown in Figure [Ml Here, the data are divided into three 
Mi luminosity intervals and compared with the data for 
the same ranges in the six main program galaxies that 
have the largest GC populations. In each case the "back- 
ground" histograms have been multiplied up by a factor 
of 6 to match the total area coverage. 

The presence of GCs in the field sample can be seen 
most noticeably in the upper two panels (A// < —9.4), 
where low peaks appear in the shaded histograms at 
{B — I)o ~ 1.6 and 2.0, near the normal positions of 
the blue and red sequences. Otherwise, the field popula- 
tions are consistent with a low and roughly flat, uniform 
distribution in color across our main color range of inter- 
est {B — J)o = 1.3 — 2.3 where the GC population lies. 
Subtraction of a uniform pedestal of only a few percent 
amplitude does not change the identification of the GC 
sequence peak colors in the bimodal fits, and has the dis- 
advantage of simply adding more noise to the data. The 
histograms have therefore not been explicitly corrected 
for this small and near-uniform background. 

6.3. Characterizing the MMR 

The next step in evaluating the color/metallicity dis- 
tributions is to characterize their detailed structure. In 
Fig. 1211 the solid lines come from the fits to the en- 
tire database and thus are weighted towards the galaxies 
with the biggest GC populations (in decreasing order of 
importance, they are NGC 4696, 7626, 3258, 3268, 1407, 
and 3348; the first two make up half the total and the last 
two only 17%). By comparison, the individual points in 
the figure give equal weight to all galaxies. Nevertheless, 
both approaches form a consistent pattern. The blue se- 
quences in the galaxies show a trend to become redder 
at brighter magnitude and more strongly so along the 
upper half of their run. No equally strong trend appears 
along the red sequences. 

Before attempting to parametrize either sequence, we 
first note that the most obvious transition region ap- 
pears to be in the range — 10 > Mj > — 11, below which 
no significant MMR appears along either sequence, and 
above which the blue sequence is much closer to the red 
sequence: 

• Along the blue sequence, the mean color for the 
three brightest bins in Table [6l for Mj < —10.4 
is {B - I)o = 1.689 ± 0.035. For the five fainter 
bins Mi > -10.4, we have {B - I)o = 1.573 ± 
0.009. The color difference between the brighter 
blue clusters and the fainter ones is thus A(_B — 
/) = 0.116 ± 0.036, which is significant at the 3.2- 
sigma level. A standard Kolmogorov-Smirnov test 
between these two histograms shows that they are 
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significantly different at much higlier than the 99% 
confidence level. The largest deviation between the 
two cumulative distributions comes at {B — I)o = 
1.6, precisely in the middle of the blue sequence. 

• Along the red sequence, {B — I)q = 1.946 ± 0.026 
for the bright interval Mj < -10.4, while {B - 
I)q = 1. 985 ±0.009 for M/ > -10.4. The difference 
A{B - /) = -0.039 ± 0.028 is significant at only 
the 1.4-sigma level. 

Another way to illustrate this luminosity-based trend 
is shown in Figure [25l and return s to the ver y first papers 
on the subject, bv lOstrov et al.l (l998) and D irsch et al.l 
((2003,). These authors used GC photometry in the For- 
nax giant NGC 1399 to suggest that the brightest clus- 
ters made up a broad and perhaps unimodal color dis- 
tribution, seemingly counter to the bimodality paradigm 
that had already become well established. In Figure \7E[ 
the color distributions over these two biggest magnitude 
intervals are shown for our six BCGs, along with the 
double-Gaussian deconvolutions for each one. The quan- 
titative RMIX fits show that (a) for the brighter interval 
Mj < —10.4 where the two modes are more strongly 
overlapped, a unimodal fit is nevertheless strongly re- 
jected in favor of two components; (b) only two modes 
are necessary for an excellent fit to any part of the mag- 
nitude range; that is, the GC population in these galaxies 
is bimodal rather than multimodal. 

As a demonstration that these fits are not dependent 
on the particular magnitude being used. Figure [57] shows 
the same comparison where now the sample is sorted by 
absolute B magnitude instead of /. The redder com- 
ponent is less prominent here and less certainly fit be- 
cause the red-sequence GCs are much fainter in B than 
in /, but the same two modes are present with very 
much the same color difference as in the M/— based 
fits. For the bright interval Mg < —9.0, we have 
/xi = 1.650±0.022, ^2 = 1.886±0.030, while for the faint 
interval Mb > -9.0, we have /xi = 1.556 ± 0.004, ^2 = 
1.971 ± 0.006. The fits in Mi are to be preferred, first 
because the /—band is much more represe ntative of the 
true b olometric luminosities of the clusters (|Mieske et al.l 
l2006f) . and second because (in the current dataset) the 
/—band photometry is slightly deeper and of higher pre- 
cision. 

Perhaps the cleanest version of the current GC data is 
presented in Figure [28l Here, the 7831 datapoints from 
all six galaxies, for the objects that have rt > 1.5 pc, 
are shown for the purpose of selecting a "best" dataset 
that minimizes field contamination, as discussed above. 
The MMR along each of the two components is grad- 
ual enough that it can be conveniently represented by 
simple polynomial fits. Useful interpolation equations in 
{Mi,B~ I) are 

(B-/)o = (2.543±0.513) + (0.237±0.103)M/ + (0.0143±0.' 

(5) 

for the blue sequence, and 

{B - /)o = (2.142 ± 0.089) + (0.017 ± 0.009)A//. (6) 

for the red sequence. These are derived from the mean 
points in Table [H] and are shown in Fig. [55] as the solid 
blue lines. The colors can be converted to metallic- 
ity through the standard transformation from the Milky 



Way clusters ([Harris et al.ll2006f) . 

[B - I)o = 2.158 + 0.375 [Fe/H] . (7) 

However, these curves should not be used to extrapolate 
to fainter levels Mj > — 8. 

The red sequence has a barely significant slope, and is 
scarcely distinguishable from a vertical line at constant 
color at {B — I)o = 1.97. In fact, the brightest two points 
along the red sequence are the least well determined in 
the whole dataset because of their relatively low pop- 
ulation. Very minor differences in the adopted RMIX 
constraints (or lack of them) can shift the mean colors of 
these two bins by a few hundredths of a magnitude.'^ 

An alternate form of the interpolation equations can 
be estimated more crudely by direct polynomial fits to 
the individual datapoints, bypassing the RMIX bimodal 
solutions entirely. To do this, I define the blue GCs to 
be the 4110 objects that fall within the color range 1.3 < 
{B - 1)0 < 1.8, and the red GCs the 3625 objects that 
faU within 1.8 < (/? - /)o < 2.4. Then the polynomial 
fits become 

(B - /)o = 2.186 + 0.152A// + 0.0094^/^^ (g) 

for the blue sequence, and 

(B - /)o = 2.183 -I- O.OI8A/7 (9) 

for the red sequence. These are shown as the dotted lines 
in Fig. [28] The uncertainties on the fitted coefficients are 
the same as for the previous pair of equations. This sec- 
ond pair of interpolation equations generally agrees well 
with the first pair, but differs most noticeably at the 
highest-luminosity part of the diagram where the data 
are sparsest. The key potential problem with the second 
approach is that it ignores the effects of mutual contam- 
ination of the two sequences on each other: that is, some 
points belonging to the red sequence scatter blueward of 
the dividing line at {B — I)o = 1.8, and conversely some 
blue-sequence points fall redward of the same line. If the 
relative numbers of objects and the intrinsic widths of 
each sequence are not equal (which is certainly the case 
here), then the amounts of mutual contamination will 
not be the same and may bias the solutions for the equa- 
tions. For this reason, I recommend that the first pair of 
equations (5,6) (the ones starting from the mean RMIX 
points) should be preferred because they more accurately 
pick out the centers of each sequence. 

6.4. Discussion and Comparisons 

To summarize this part of the analysis, the very most 
basic conclusion is that the mean color of the blue- 
sequence clusters begins to change, with high statisti- 
cal significance, at around the luminosity level Mj ~ 
— 10.5. For a typ ical mass-to-light ratio {M/L)j = 2 
(iMcLaughlinl [2000') . this level translates to a GC mass 
'MMjIO^Mq. Below this level, the present data give no 
compelling evidence that there is a significant correlation 
between luminosity and color. Above this level, the blue 

^ See also IMleske et al.l II2006I ). who demonstrate that the 
bimodal-Gaussian fitting procedure can end up producing a slight 
artificially positive slope to the red sequence if the two sequences 
are partially overlapped and not equally populated; see their sec- 
tion 3.2. The fitted slope deduced here is consistent with their 
results. 
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GCs become significantly redder in the mean, causing 
the blue and red modes to overlap more strongly. 

At the same time, this transition point does not seem 
to be abrupt. Figs. [20l and [21] suggest that the blue se- 
quence makes a gradual changeover from a near-vertical 
sequence at fainter magnitudes, to a sloped one at 
brighter levels. The polynomial fits stated above reflect 
this changeover. For the range M/ < —9.5 crossing the 
transition region, if we simply approximate the trend 
by a linear fit, the mean slope is A(i? — /)/A(M/) = 
-0.072 ± 0.018. Since A{B - /)/A([Fe/H]) = 0.375 
(|Harris et al.l l20d6h . this translates to a heavy-element 
abundance scaling Z L^^**"'^^. Extrapolating up- 
ward to the very most massive GCs known, this corre- 
lation suggests that the blue-sequence GCs at lO^M© 
should have an average (Fe/H) ~ —1.0 instead of the 
"normal" mean [Fe/H] ~ —1.5 which the vast majority 
of the lower-mass ones have. 

An "anomalous" feature in the bright-end distribution 
appears in the Mj range around —10.5 to —11, where 
in Fig. [m three of the red-sequence points stand notice- 
ably off the main distribution. These three points are at 
{B — I)o ~ 1.8 and are all from the two galaxies NGC 
3258 and 3268 (see Fig. Referring back to the raw 
CMDs fFigs. [T^ and fTH]) we see that this is the magni- 
tude range where the overlap and blurring between the 
blue and red sequences becomes quite noticeable. The 
heavier overlap decreases the stability of the bimodal fit- 
ting and increases the internal uncertainties in all of the 
fitted parameters (/ii,2, o-i,2,Pi/p2)- My provisional in- 
terpretation is therefore that these anomalous points are 
primarily the result of small-sample statistics; the fact 
that these intermediate-color points do not show up in 
the bimodal fits to the combined sample of all six galax- 
ies with its higher statistical weight (Table [6|) provides 
some support for that conclusion. Nevertheless, the pos- 
sibility should be kept in mind that our basic assumption 
of a two-sequence model - the assumption that is built 
in to almost every previous discussion of GC metallicity 
distributions - may not capture the entire nature of the 
color/metallicity distribution for these extremely high- 
mass clusters. 

Because the photometry of the program objects is 
specifically corrected for their individual effective radii, 
any concerns that the basic existence of the MMR is 
merely an artifact of the measure ment process are re - 
moved. This study, and the one bv iMieske et"aLl (|2006D 
for the Virgo galaxies, now provide the major HST-based 
datasets in which the GC photometry is accurately size- 
corrected. 

Extensive comparisons of the different forms of MMR 
published in earlier papers are beyond the aims of the 
present work, though more su ch intercomparisons are 
done by ICockcroft et al.l ()2009D . They combine the re- 
sults from most of these previous studies by convert- 
ing the several different photometric color systems that 
have been used for different galaxies into a homoge- 
neous metallicity scale, and correlate the results for the 
MMR slopes with galaxy luminosity and GCS specific 
frequency. 

A still newer study by [Waters et al.l (|2009| ) of the sin- 
gle galaxy M87 is of particular interest for the unusual 
photometric characteristics of the data. They use ex- 



tremely deep HST/ACS imaging in V and / covering the 
inner 10 kpc (1.5i?e//) of M87 to construct aperture- 
corrected photometry of about 2000 GCs, deep enough 
to cover almost the entire GCLF. Of these, ~ 1000 are 
brighter than the GCLF turnover point (compared with 
more than 8000 from the six galaxies in this paper) and 
about 130 are brighter than Mj = —10.5 (compared with 
over 750 from the present study) . They find the normal, 
first-order bimodality of the MDF although with a large 
internal spread in the {V — I) colors, but also claim that 
that the blue GC sequence has no significant MMR slope. 
This conclusion relies on a forced linear-fit model over the 
whole magnitude range in (J, V — I) and is thus domi- 
nated by the 94% of the sample that is fainter than the 
actual MMR "turn-on" point at Mj ~ -10.5. In fact, 
their data do not rule out the type of MMR that has been 
found here: that is, a near-vertical blue sequence on the 
fainter part of the magnitude range, moving upward to 
a stronger slope towards higher metallicity with a tran- 
sition region above a million Solar masses. The data in 
the present study draw from a statistical sample more 
than 6 times larger, and in addition, use a photometric 
index {B — I) that is twice as sensitive to GC metallicity 
as is (V -I). 

Of the previous studi es, the most closely comparable 
one is probably Micsk e et al.l (|200 6) . who analyze GC 
samples from 79 Virgo ellipticals and perform bimodal- 
Gaussian fits to the GC color distributions with tech- 
niques very similar to the ones used here. Their database 
has similar size (a total of a few thousand GCs brighter 
than the GCLF turnover luminosity), uses a color index 
(g' — z') with similar metallicity sensitivity, and has a low 
level of field contamination. They divide the galaxies into 
four different luminosity groups (divided roughly speak- 
ing into supergiants, giants, normal ellipticals, dwarfs) 
and then discuss the trends of mean GC color versus 
magnitude in each group. The most important differ- 
ences between their analysis and the present one are that 
(a) the Virgo data extend about 1.5 mag fainter in GC 
luminosity, below the GCLF turnover; and (b) they ap- 
ply the forced linear-fit assumption to the mean points 
from the bimodal fits, which I have argued against above. 
The blue GC sequence is found to have a significant slope 
(MMR) for the three bins containing the large ellipticals, 
but no significant slope is found for the GCs in the dwarf 
ellipticals. In all cases the red GC sequence is indistin- 
guishable from vertical. Referring particularly to their 
Figure 2 where the mean points are plotted, I note that 
the significance of the blue MMR slope will be very much 
reduced if the brightest two points (M^ < —10, close 
to the IO^Mq level) are excluded and only the fainter 
bins are used. The strongest case in their data for a 
nonzero MMR slope is from the combined GC sample 
for the three brightest giants, M87, M49, and M60 (their 
Fig. 2a), though even there the blue sequence becomes 
essentially vertical for > — 9. 

Additional tests for any biasses and internal uncer- 
tainties in the fitting procedures can be carried out 
through simulated CMDs (see, e.g. iMieske et all 120061 : 
IWaters et al.l l2Q09l for recent examples). Figure [29] 
shows one illustration of a model CMD tuned to the 
present {B — I) data, in which the total numbers of 
GCs and contaminating field objects are closely similar 
to the real data for the six BCGs combined. Here, the 
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simulated blue sequence is purely vertical with a mean 
color {B — I) = 1.56 and Gaussian standard deviation 
as-i — 0.09, and follows a Gaussian luminosity function 
with (Mj) = —8.3, ai = 1.45. The red sequence is purely 
vertical with {B - I) = 2.00, as-i = 0.15, and its LF 
follows (Ml) = —8.3, (7j = 1.25. Field stars are added 
following a uniform distribution in color and a power- 
law increase in magnitude. All the artificial data are 
broadened in color by measurement scatter that increases 
with magnitude, and lastly have photometric complete- 
ness functions in B and / applied to them that mimic 
the real data. The simulated CMD in Figl29] includes 
17500 GCs, 55% of which are on the red sequence and 
45% on the blue sequence. Lastly, the histograms above 
each CMD show that the intrinsic distributions in (B — I) 
for the luminosity range below the onset of the MMR are 
similar. The simulated CMD clearly fails to reproduce 
the structure at the high-luminosity end (M/ < —10) 
that we see in the real data. The most noticeable effect 
is that the real data do not have the very metal-poor and 
very luminous GCs that would lie on the direct upward 
extension of a vertical blue sequence. If these were actu- 
ally present in the real galaxies, they would be more eas- 
ily found than GCs anywhere else in the color-magnitude 
diagram. RMIX bimodal fits to the simulated data show 
no significant trend for a spurious MMR to appear, or 
for significant biasses in the deduced sequence colors. 

6.5. Implications for Theoretical Models 

The currently published analyses that draw from large, 
deep, and metallicity-sensitive photometry point to an 
emerging consensus around the view that the MMR has a 
nonlinear form: the blue sequence is nearly vertical at low 
GC luminosities, but continuously changes to a stronger 
increase of metallicity with cluster mass as we go upward 
along the sequence. On the red sequence, however, the 
current data are consistent with the interpretation that 
it stays nearly vertical (metallicity is uncorrelated with 
luminosity) . 

These observations, if correct, already have implica- 
tions for the different theoretical mo dels based on the 
centra l concept of self-enrichment. IStrader fc SmithI 
(|2008l ) construct a model built on the assumption of 
a continuous linear form of the MMR, which as noted 
above corresponds to a scaling Z ~ over the en- 
tire GC mass range. In their model this form turns out 
to require that the star formation efficiency within 
the protocluster cloud also scales continuously with clus- 
ter mass in the same way. This particular form of a 
self-enrichment model is therefore ruled out by the ob- 
servations, which show that the blue sequence is nearly 
vertical (metallicity is uncorrelated with mass) for M < 
1O^M0, and t hat more genera l ly the MMR is nonlinear. 

In contrast. iBailin fc Harris! (|2009[ ) assume — const 
with GC mass (specifically = 0.3, though the pre- 
cise value is not important). The important quantity 
that changes with cluster mass is instead the fraction of 
SNe-ejected gas held back in the protocluster, i.e. the 
heavy-element retention efficiency fz which drives the 
enrichment of the lower-mass stars in the protocluster. 
The retention fz is a strong nonlinear function of M 
(that is, the depth of the protocluster potential well). 
This approach leads to the prediction that the blue GC 
sequence should show a fairly distinct visible "onset" of 



the MMR at around one million Solar masses. The same 
sort of self-enrichment onset should also appear on the 
red GC sequence, though it would only become obvious 
at still higher mass because the mean pre-enrichment of 
the red GCs (that is, their baseline mean metallicity) is 
an order of magnitude higher. This theoretical form of 
the MMR is more in agreement with the observations. 

Both these model approaches are admittedly ideal- 
ized, and a synthesis of their ideas might lead to an im- 
proved match to the data. The Bailin/Harris model pre- 
dicts a MMR onset that is probably too abrupt (a rapid 
changeover to a sloped MMR appears just at the point 
where the amount of internal self-enrichment becomes 
larger than the baseline pre-enrichment metallicity) . The 
data here (see particularly Fig. [2T|) favor a smoother and 
more gradual change in slope over a wide mass range. 
Clearly the enriched-gas retention efficiency fz should 
increase with protocluster mass, but the star formation 
efficiency may also change systematically with increasing 
mass. These models also treat the protocluster in isola- 
tion. If it is placed within a surrounding lower-density 
gaseous region (giant molecular cloud or dwarf galaxy), 
then the addition of some external pressure confinement 
could increase the effective fz and extend the visible on- 
set of the MMR more smoothly to lower masses (see the 
model papers for additional discussion). 

6.6. The MMR and Cluster Size 

A final stage of the MMR discussion in this paper is to 
ask whether the details of the mass/metallicity correla- 
tion might depend on other external factors. Two factors 
that can be investigated here are the characteristic size 
of the GC and its location in the galaxy, both of which 
may affect the efficiency of self-enrichment. 

If the MMR itself is a second- order feature of the 
GC metallicity distribution, requiring 10'^ clusters to 
clearly reveal it, then any correlations of the MMR with 
other parameters are still more subtle third- order fea- 
tures that we might expect to require samples of ~ 10"* 
clusters or more. The database in the present paper is 
therefore just barely large enough to begin this kind of 
search. Indeed, the results to be described below are 
suggestive of certain trends but in the end must be con- 
sidered as no more than tantalizing hints. 

The slope of the MMR should in principle depend on 
the depth of the cluster potential well, and its effective 
radius helps determine that depth. That is, for a given 
cluster mass we might expect the self-enrichment process 
to be less effective for more extended GCs, which would 
have lower escape energies and thus be less able to hold 
back their early SNe ejecta. If this idea is correct, then 
along the blue GC sequence the clusters with the biggest 
radii should be systematically bluer (more metal-poor) 
than those with small radii. 

A first look for any such effect is shown in Figure [30l 
where the total sample is subdivided by r^. The only 
pre-selection applied here is to restrict the list to those 
objects with ISHAPE S/N > 20 to ensure some minimal 
quality in the cluster size estimate. 

• As discussed above, the subsample with r/i < 2 pc 
is likely to have the highest proportion of field con- 
taminants of various kinds. Nevertheless, many of 
the objects in this range are likely to be GCs: the 
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bimodal sequences remain well defined and their lo- 
cations track the mean loci defined from the whole 
sample (dashed lines in the Figure). 

• The objects in the range 2 < r?i < 4 pc make 
up the classic, "normal" globular clusters whose 
sizes fall comfortably into the baseline pattern set 
by the Milky Way. This part of the sample is 
likely to be relatively free from contamination, 
and defines clear bimodal sequences. They accu- 
rately follow the mean lines for the entire sam- 
ple as well. Interestingly, there are relatively more 
/ii(7/i-luminosity clusters (M/ < —11) in this group 
than in the previous one, but this is largely due to 
the radius/luminosity correlation discussed earlier 

(Fig.m- 

• The biggest objects are in the third panel of the 
figure. These fall in the size range 4 < r/i < 10 pc 
and would include the most extended known GCs 
with luminosities higher than the GCLF turnover 
point, except for a few still bigger UCD-like ob- 
jects. One feature to note is that, at all luminosi- 
ties, there are relatively fewer red-sequence clus- 
ters in this size range than at the lower sizes; this 
is the result of the size/metallicity correlation dis- 
cussed earlier (Fig.[T8|). The details of the blue GC 
sequence are potentially more interesting, because 
if the predictions of the self-enrichment model are 
correct, we might expect to begin seeing its effects 
here. Indeed, these large clusters fall steadily about 
0.1 mag to the blue of the mean line for the whole 
sample, corresponding to A[Fe/II] ~ —0.25 dex in 
metallicity. The self-enrichment model is likely not 
to be the only possible explanation of this offset, 
but it is a tantalizing piece of evidence that at very 
least is consistent with the model. 

A Kolmogorov-Smirnov test was used to compare the 
color distributions for the three groups of clusters shown 
in Fig.[30l restricted to the upper luminosity range Mj < 
— 10.0 where the differences appear to be most noticeable. 
The third group (rh > 4 pc) is systematically bluer than 
the first group {r^ < 2 pc) with very high significance, 
> 99%. The middle group {rh = 2—4 pc) is different from 
either of the other two at the > 90% level, suggestive but 
not conclusive. 

In the third panel of Fig. [SOF c). the very high lumi- 
nosity range —11.2 > Mj > —11.8 also shows a dozen 
objects that fall between the two main sequences. Two 
are from the NGC 3268 field, seven from NGC 7626, and 
three from NGC 4696, aU have S/N > 200, and all have 
Th values from 4 to 7 pc. Careful visual inspection on the 
original images shows that one of these has a comparably 
bright neighbor 0.45" away, not close enough to invali- 
date the profile fits. All the others are isolated objects 
with small, GC-like ellipticities and no morphological pe- 
culiarities. There is no reason to reject them from the 
GC sample. 

6.7. Trends versus Galactocentric Distance 

iMieske et al.l (|2006l l , in their study of the Virgo galaxy 
GCs, suggested that the MMR may become less notice- 
able at larger galactocentric distance. Such a trend might 



be the mark of different enrichment or formation his- 
tories at different locations in the parent galaxy's halo. 
However, their data were restricted to the inner radial 
region Rgc < 1.5i?e// (for the bigger galaxies). The 
present sample covers a much bigger range in Rgc and 
thus should be able to provide a more sensitive test of 
any such trend. Figure [3T] shows the color-magnitude di- 
agrams in three successive radial zones divided by galac- 
tocentric distance: the inner region Rgc < Reff , a middle 
zone Reff < Rgc < '2Reff, and the broad outer region 
^Reff < Rgc ^ ?>Reff- These radial subdivisions leave 
about the same numbers of objects in each zone, so the 
mean red and blue sequence lines that were defined pre- 
viously (the heavy dashed lines in the figure) represent 
essentially a straight average over all radial zones. 

The correlation of MMR with Rgc also provides an 
implicit secondary test that the GC photometry is not 
systematically affected by the background light gradient 
of the galaxy. For Rgc > Reff, the gradient becomes 
low and any residual effects on the measurements will be 
small. 

To search for clues that the MMR slope might depend 
on location, the RMIX bimodal fits described above were 
repeated separately for each zone. The results are shown 
in Figure [32l Here, because the sample sizes in each lu- 
minosity bin are smaller, it was not possible to obtain 
fully unconstrained fits to the same level of confidence 
for all five parameters (/ii, /i2, txi, cr2,pi/p2) in the high- 
luminosity, sparsely populated bins that are the key to 
the MMR slope. However, because we are looking es- 
pecially for differences in the mean colors of the blue 
sequence from one zone to the next, the RMIX fits were 
deliberately constrained so that the mean color and stan- 
dard deviation of the red sequence were fixed, since it is 
nearly vertical to begin with. The standard deviation is 
fixed at a2{red) — 0.18 (see TableE]), and the mean color 
fixed at fj,2{red) = 1.994 for the inner zone, 1.979 for the 
middle zone, and 1.949 for the outer zone. These values 
were adopted to account for the shallow radial metallic- 
ity gradient in the system (discussed in the next section) . 
In addition, to smooth the mean sequences numerically, 
the RMIX fits were done in ~ 0.5— mag wide luminosity 
bins but stepped by AM/ = 0.25 mag between bins. 

In Fig. [32l from bottom to top the same gradual 
changeover as noted previously is apparent, from a nearly 
vertical sequence at Mj < —10 to a progressively steeper 
MMR at higher luminosities. The internal errorbars on 
each mean point also grow dramatically at higher lumi- 
nosity as the sample size in each bin decreases; and for 
M] ^ —11, the bimodal fits become rather unstable as 
the blue sequence overlaps more and more heavily with 
the red sequence and even the proportions {pi , P2) in each 
component become quite uncertain. There is no convinc- 
ing difference between the mean lines for the inner (solid 
line) and middle (dashed line) zones. A more noticeable 
difference appears for the outer zone, which stays offset 
to the blue by ~ 0.1 mag relative to the inner zones. 
Taken at face value, this means that the MMR decreases 
in amplitude further out in the halo. 

Another test that avoids some of the problems of bin- 
ning the data, is to look directly at the color distribu- 
tion in {B — 1)0 for the blue-sequence GCs brighter than 
M] = —10.5. The fiducial dividing line between the 
red and blue sequences is {B — /)o ~ 1.8: the inner 
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zone has 97 objects bluer than that boundary, the mid- 
dle zone 139, and the outer zone 129. Figure [33] shows 
the cumulative color distributions for these bright, blue 
GCs. The existence of a global, weak radial gradient in 
mean metallicity (see next section) first requires adding 
A{B — I) = 0.02 mag to the middle zone and 0.03 mag 
to the outer zone before comparison with the inner zone. 
The color distributions becomes steadily more weighted 
to the bluer side at larger R/R^ff, particularly in the 
color range [B — I)o ~ 1.6 — 1.7 that is crucial for de- 
termining the MMR. However, a normal Kolmogorov- 
Smirnov two-sample test shows that none of the zones 
are different from any of the others at higher than 70% 
significance. This is not sufficient to demonstrate any 
radial change in the MMR convincingly and is left only 
as another tantalizing hint. 

A tentative conclusion from these tests is that (a) the 
MMR exists at all the galactocentric distance zones that 
the present data are capable of exploring; but that (b) 
there is very weak evidence that the MMR becomes grad- 
ually less noticeable at larger galactocentric distance. 
The directio n of t his effect is the same as found by 
iMieske et all ()2006l ) (see their Table 2 and section 4.2 
discussion). Bigger samples of high-luminosity clusters 
from other galaxies will be needed to probe this effect 
c onclusively. 

iBailin fc Harris! ()2009l ) suggest two effects that might 
end up producing a radial gradient in the MMR. One 
of these would be relevant in large galaxies that hosted 
active galactic nuclei (AGN's) at an early stage in their 
evolution: the AGN would dramatically increase the ex- 
ternal ultraviolet radiation field and, possibly, inhibit the 
self-enrichment process within a protoglobular cluster of 
a given mass. Taken at face value, this model would pre- 
dict that the MMR should be stronger at larger galacto- 
centric distance farther from the AGN, which is opposite 
to what we observe. However, there is no guarantee that 
the timing of GC formation would automatically overlap 
with the appearance of the central AGN. The UV radi- 
ation field might even produce the opposite gradient if, 
for example, the inner GCs had already finished forming 
by the time the AGN reached its peak luminosity but the 
outer GCs had not. 

The second possible effect is related to the difference 
between the initial masses of the GCs (which determine 
their degree of self-enrichment) and their present-day 
masses, which are smaller because of mass loss due to 
early SNe ejecta, stellar winds, and the slower processes 
of tidal stripping and evaporation of stars over 10^" y. 
The inner GCs should have been subjected to a higher 
degree of dynamical mass loss due to tidal stripping than 
the outer GCs over their lifetimes, so an inner-halo GC 
that has the same mass now as an outer-halo one would 
have started out at higher mass and thus have been able 
to self-enrich more. Testing this idea more quantitatively 
and finding out specifically how it affects the final MMR 
will require knowing, at a minimum, the average amount 
of GC mass loss from the protocluster stage onward, as 
a function of both MGc{i™t) and Rgc- This direction 
may prove to be a rich area for detailed modelling. 

Another GC phenomenon that may be linked to these 
arguments is the long-standing puzzle that the GC lumi- 
nosity (or mass) function is nearly independent of galac- 
tocentric distance r. If the initial GCMF were nearly the 



same throughout the early galactic halo, then the strong 
dependence of dynamical evolution time on r should 
build in a present-day GCMF whose low-mass end and 
"turnover point" depend on r. A recent discussion by 
iMcLaughlin fc Fall! (2008) shows that the observed shape 
of the GCMF can be explained if its initial shape depends 
appropriately on the characteristic cluster density ph, be- 
cause the mass loss rates increase with p. The fact that 
the outer clusters have lower average densities would also 
go along with less self-enrichment and a weaker MMR at 
larger distances from the galactic center. 

7. METALLICITY GRADIENTS 

Careful inspection of the mean positions of the near- 
vertical sections of both blue and red sequence lines 
shows that at larger radius, the mean color of each se- 
quence shifts slightly to the blue compared with the loci 
at smaller radu. That is, both types of GCs appear 
to show a metallicity gradient. It has long been real- 
ized that the GCS in a normal large galaxy taken as 
a whole (that is, combining all types of clusters) shows 
a negative metallicity gradient, but the majority of the 
effect is created simply by the decreasing ratio of red 
versus blue GCs w i th increasing galactocentric distance 
(iGeisler et al."1996': 'Harris et al. 1998; Dirsch eLaL||2003|; 
iRhode fc Z cpl 2004; Harris et al. 2006). That is, the 
steeper spatial distribution of the red, metal-rich clusters 
produces a population gradient, whic h has the byproduct 
of an overall metallicity gradient. In [Harris et al.l (j2006f ) 
we discussed the population gradients in the BCGs stud- 
ied here. By contrast, much less published material exists 
that is directed at finding any true metallicity gradient 
within each of the two types of clusters, and indeed any 
such gradients are subtle. 

A large body of literature exists for abundance gra- 
dients in the integrated light of large galaxies, mea- 
sured from line-strength indices. These studies are usu- 
ally based on data from the high-surface-b rightness re- 
gion within IReff of the gala.xy ce i iter (seelFisher et alj 
iggs"; 'Kobavashi fc Arimotdjllli iMehlert et al.l l2003t 
Sanchcz-Blanquez ct al. 200Q, for a sampling of such 
studies, along with extensive guides to the literature). 
The analyses typically show that heavy-element abun- 
dance within this bulge light scales with projected galac- 
tocentric distance as anywhere from Z ^ R~'~'-^^ to 
i?"*^'^, favoring interpretations in which the central re- 
gions of these galaxies formed by a combination of dis- 
sipative collapse (favoring steep gradients) and subse- 
quent mergers (favoring flatter gradients). The results 
are rather strongly model-dependent in the sense that 
each line index must be transformed to some tracer of 
heavy-element abundance. 

Globular clusters can provide valuable extensions to 
these integrated-light measurements. GCs are easy to 
measure at radii well beyond Reff, and transformation 
of their colors to metallicities is based directly on calibra- 
tion from the Milky Way clusters. These transformations 
need only the reasonable assumption that the GCs are 
"old" (> 10 Gyr) and thus in the regime where the colors 
are quite insensitive to minor age differences. The only 
practical issue is that a large statistical sample of GCs, 
over a large radial range, is needed to measure the metal- 
licity gradient accurately. A large radial range, and large 
sample size, are precisely the advantages of the database 
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in the present study. 

Figure [34] shows the intrinsic colors {B — I)o of the GCs 
in all six galaxies plotted directly against galactocentric 
distance, again normalized to Ref f for each one. For this 
purpose the sample is restricted to the brighter objects 
{Mj < —9) with positive size measurements {rh > 0), 
to minimize the effects of field contamination. As above, 
the dividing line between blue and red GCs is defined 
as {B — I)q = 1.80 and we then solve directly for the 
best-fitting relation for mean color versus radius within 
each group, in the conventional form {B — I)q = a + 
b\og{R/R,ff). 

Table [7| lists the solutions for the individual galaxies, 
along with the solutions for the combined sample. The 
slopes b in all cases are shallow, and essentially equal for 
both the red and blue components. In addition, they are 
not significantly different between individual galaxies, so 
to gain statistical weight, we can combine all of them into 
a single solution. The combined sample yields a nonzero 
signal for the slopes at the > 4.5 — a level for both types 
of clusters, and therefore the existence of a metallicity 
gradient appears to be convincingly real. Given the con- 
version ratio A{B-I) = 0.375A[Fe/H] (see above), both 
of the gradients convert to a heavy-element abundance 
scaling Zgc R-°-^°±°-°^\ 

The slope AlogZ/Alogi? = —0.10 is shallower than is 
conventionally found in the integrated-light studies (see 
the references cited above), but it is also dominated by 
more distant regions of the halo, R > Reff- Few other 
GCS studies exist for whic h direct comparisons can be 
made. Geisler ^t al.l (|199 6f), in their wide- field study of 
the GCS around the Virgo giant NGC 4472, found slopes 
of Z ~ i?-o.i5±o.o3 foj. ^Yie blue GCs and Z - i?-o i2±o.o6 

for the blue GCs, out to a limiting radius of ~ BReff. 
These are clearly similar to the results here. Geisler et 
al. also point out that the gradient for the red GCs 
closely matches the color gradient for the integrated halo 
light of the host galaxy; and dire ct evidence from halo- 
star studies in giant e llipticals (e.g. [Harris fc Harrisl[200l 
iReikuba et al.l [2005! ) shows that the field-star popula- 
tion is dominated by the metal-richer stars whose MDF 
matc h es that of the red GCs (see also lForte et al.|[2005l 
|2007( ). iForte et all ()2001[ ) determined the gradient for the 
blue GCs in two Fornax ellipticals, NGC 1399 and 1427, 
finding Z{blue) ~ R~^-^ for both systems. Their data 
extended to 12 ^e ff for NGC 13 99 and lOReff for NGC 
1427. Recently, iLee et al.l ([2008) have found shallow but 
significant (C — Ti) color gradients in both the blue and 
red GCs in the Virgo giant M6 that are quite simil ar to 
those for M 49. Lastly, in M31. lPerrett et al.l ()2002f ) and 
iFan et aTl (2008) find that the metal-poor GCs exhibit a 
metallicity gradient out to 20 kpc (or about AR^ff) but 
suggest that the metal-rich GCs show no gradient. 

Perhaps the ea r hest m odern investigation of this ques- 
tion is by IZinnI (|1985l ) for the Milky Way itself. He 
found that the inner, metal-poor GCs were significantly 
more enriched than the outer-halo ones. Using contem- 
porary metallicity data for 105 metal-poor clusters with 
[F e/H] < -0.95 (with data taken from the 2003 edition 
of lHarrisl[l9M ). I find a best-fit solution ( Fe/H ) = 
(-1.446 ±0.063)- (0.165 ±0.062)logrgc where rgc is the 
true three-dimensional distance from the Galactic center. 
This GC sample extends out to more than 100 kpc from 
the Galactic center. A solution in which the same data 



are artificially projected onto two dimensions (making 
them directly comparable with the other galaxies above) 
gives ( Fe/H ) = (-1.490±0.047) - (0.146±0.051)logi?gc. 
In this case the 3D slope is already so shallow that pro- 
jection into two dimensions does not decrease it much. 

The available data sketch out a consistent pattern 
whereby both the metal-poor and metal-rich globular 
cluster subsystems have shallow but real radial metallic- 
ity gradients, in the range 

Z ^ i?-(o i-o.2)^ Particularly 
in the large spiral galaxies, the gradients for the metal- 
rich clusters are harder to determine because they are 
less numerous and spatially more centrally concentrated, 
but for the metal-poor clusters the Milky Way scaling 
of Z ~ /j-o.i5±o.05 jg typical of the mean for aU the 
galaxies. The most internally precise value for any of the 
systems is the one determined here (6 — —0.10 ± 0.022) 
for the six BCGs. 

Interpretation of these metallicity gradients in terms 
of the past history of the parent galaxy is an intriguing 
and incompletely explored area. Major mergers can pro- 
duce extensive spatial mixing and reduce any pre-existing 
steep metallicity gradient that was the result of dissipa- 
tive early star formation in the progenitors. Gas-rich 
mergers will produce a new population of stars and clus- 
ters especially in the inner regions and may thus help 
rebuild the gradient, at least for the metal-rich com- 
ponent. However, two intriguing features of the obser- 
vational evidence so far are that (a) the gradient slope 
b ~ —(0.1 — 0.2) appears to be very much the same for 
both giant elliptical and spiral galaxies, and (b) it is also 
the same for the metal-rich and metal-poor clusters, even 
though they may have experienced very different histo- 
ries as well as formation times. More extensive N-body 
simulations of mergers specially tuned to trace metallic- 
ity gradients well outward into the halo will be needed to 
constrain the range of possibilities and help decide how 
much of the gradients we observe now are relics of the 
primordial conditions of formation. 

8. SUMMARY 

In this study the globular cluster populations in six 
giant elliptical galaxies have been investigated through 
their distributions in size, metallicity, and luminosity. 
A total of about 15000 objects, of which about ^ 8000 
are high-probability globular clusters belonging to these 
six galaxies, have been measured down to a luminosity 
Mj ~ —8, slightly fainter than the turnover point of the 
standard GC luminosity function. Particular emphasis 
is given to characterizing the relation between GC mass 
and metallicity (the MMR) that becomes most notice- 
able for the blue, metal-poor cluster population at high 
luminosity. A summary of the findings in this study can 
be made as follows: 

1. The target galaxies in this study are at typical dis- 
tances d ~ 40 Mpc and at these distances, their 
globular clusters can be barely resolved as non- 
stellar with our HST/ACS imaging data. Ana- 
lytic King-model fits convolved with the stellar 
point-spread function have been used successfully 
to determine total magnitudes, colors, and effective 
(half-light) radii for each detectable object. 

2. The new photometry of the clusters based on indi- 
vidually size-corrected total magnitudes is more ac- 
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curate than the previous data based ei ther on PSF- 
fitting or fixed-aperture photometry (jHarris et alJ 
I2OO6D . but it does not lead to substantial changes 
in the color-magnitude distribution. The GCs ac- 
curately define the classic bimodal red and blue 
components at all luminosities, with the blue com- 
ponent at {B — I)q = 1.6 and the red component 
at {B - 1)0 = 2.0. 

3. The distribution in intrinsic GC scale size in these 
giant ellipticals accurately matches that of the 
Milky Way clusters, with a peak frequency near 
r/,, ~ 2 pc and an extended tail of larger objects. 
Field contamination of the total sample consists 
mostly of very faint or starlike objects that can be 
objectively culled out of the GC data. 

4. The range of intrinsic sizes is large for all types of 
clusters, metal-poor and metal-rich alike. However, 
the blue, lower-metallicity GCs are systematically 
larger on average than the red, higher-metallicity 
ones by about 17%. This difference is statistically 
highly significant, and holds up at all GC luminosi- 
ties and all galactocentric distances. These results 
suggest that the correlation between cluster scale 
size and metallicity is intrinsic to the clusters, a 
result of their formation or later internal dynami- 
cal evolution, rather than reaction to the external 
tidal field. 

5. The median cluster sizes, of both types, also scale 
with projected galactocentric distance as ~ 
Ryc^- The data extend outward to approximately 
5Reff into the halos of the galaxies, and thus probe 
more extended regions of the halos than do most 
previous studies. These results fit in well with data 
from other galaxies indicating that clusters forming 
in shallower, less extreme tidal fields end up with 
larger characteristic sizes. 

6. The median scale size of the clusters also depends 
weakly and nonlinearly on cluster luminosity, in 
the sense that GCs more luminous than Mj < — 1 1 
become progressively bigger. The effective radii of 
the very most luminous clusters known are almost 
twice as large as the bulk of the clusters at Mj > 
— 11. Most of these luminous GCs, however, are 
smaller than UCDs of the same luminosity. 

7. The blue GC sequence exhibits a clear 
mass/metallicity relation (MMR), but it is 
not linear in color versus magnitude. Its shape 
is rather that of a gradual changeover from 
nearly vertical (for the lower-luminosity range 
Mj > —9.5) to a steeper slope at higher luminos- 
ity. The uppermost luminosity range corresponds 
to a heavy-element scaling with cluster mass of 
Z - AfO-S. 

8. The red GC sequence is nearly vertical at all lumi- 
nosities. 

9. The shape of the MMR for both the blue and red 
sequences is consistent with the basic interpreta- 
tion that it is due to cluster self-enrichment during 
their formation epoch. Though theoretical models 



of this type are still in early stages, the data rule 
out models in which the star formation efficiency 
fi, within a protocluster scales strongly with clus- 
ter mass. Instead, models in which f^, ^ const and 
the heavy-element retention efficiency fz increases 
with cluster mass produce a more realistic match 
to the observations. 

10. The slope of the MMR along the blue sequence 
does not depend strongly on cluster size r^- The 
data contain a hint that the biggest GCs {rh > 4 
pc) follow a slightly bluer sequence at the high- 
luminosity end. If so, it would be consistent with 
the idea that self-enrichment was less effective (for 
a given GC mass) in more extended clusters with 
shallower potential wells. 

11. The data from these giant elliptical galaxies also 
show that the MMR depends weakly, if at all, on 
projected galactocentric distance. If a trend exists, 
it is in the sense that the MMR has lower ampli- 
tude at larger Rgc', this may indicate that the in- 
ner GCs have lost more mass than the outer ones 
due to tidal stripping and dynamical evolution, so 
that they have a larger contrast between their ini- 
tial masses (which determine the amount of self- 
enrichment) and their present-day masses (which 
we directly observe). 

12. Both the red and blue GC subsystems exhibit shal- 
low metallicity gradients: they become system- 
atically more metal-poor at larger galactocentric 
distance. The gradients in both cases scale as 
Z Rg^-^^'^^-^^^ and, though small, are clearly 
measurable in this large dataset. These results 
agree well with the limited data available for other 
elliptical and spiral galaxies, and are likely to point 
to a combination of the amount of initial enrich- 
ment plus later merging histories. 

9. FINAL CAVEATS 

Phenomena such as the MMR, the GC size distribu- 
tions, and the metallicity gradients within the compo- 
nents of the globular cluster systems, have the potential 
to provide rich grounds for modelling and interpretation, 
but are still imperfectly understood. All the attempts to 
date to describe what is happening at the high-mass end 
of the globular cluster distribution continue to be data- 
limited. Even in the present study, only a few hundred 
GCs brighter than M/ = -10.5 (10'^ M©) are available to 
work with, and these are spread over several individual 
galaxies within which the conditions of formation could 
have differed. Breaking down that part of the sample 
further by metallicity, luminosity, or galactocentric dis- 
tance then leads us quickly into a shadowy zone where 
conclusions are tentative. 

On strictly observational grounds, the correlation be- 
tween GC color and luminosity called the MMR contin- 
ues to be a difficult effect to measure, and should still be 
considered as incompletely established. Even the data 
discussed here are not yet enough to address important 
questions such as whether or not the MMR may differ in 
amplitude from galaxy to galaxy, or if it is genuinely ab- 
sent in some galaxies. Careful attempts have been made 



Globular Cluster Systems 



19 



here to reduce or eliminate several worries in previous 
studies, including small sample size; systematic errors in 
PSF-fitting or fixed-aperture photometry due to differing 
cluster sizes; low signal-to-noise; and background light 
gradients. If the present data still hide systematic er- 
rors, they must be due to some subtle and ill-understood 
features of the image preprocessing or measurement algo- 
rithms that would preferentially affect the photometry of 
the brightest objects instead of the faintest ones, and the 
blue GCs rather than the only-slightly-redder ones. In 
the longer term, probably the best security against such 
systematics will be to carry out more GCS studies with 
different telescopes and cameras, still larger samples, and 
especially metallicity-sensitive color indices. Some com- 
parative work of this type is already in the literature, 
but there are still only a few large galaxies in which the 
MMR can be described as plainly visible. 

On theoretical grounds, the central concept of self- 
enrichment is an attractive one, but its most impor- 
tant assumption is that it requires the main star for- 
mation period in a proto-GC to last long enough for the 
most massive stars to form, evolve, and eject their en- 
riched material before most of the low-mass stars have 
finished forming themselves. The often-used arbitrary 
assumption of an "instantaneous starburst" cannot be 
correct in real star-forming regions, but for this model 



Aarseth, S.J., & Heggie, D.C. 1998, MNRAS, 297, 794 

Bailin, J., & Harris, W.E. 2009, ApJ, in press (arXiv:0901:2302) 

Barmby, P., McLaughlin, D.E., Harris, W.E., Harris, G.L.H., & 

Forbes, D.A. 2007, AJ, 133, 2764 
Baumgardt, H., Hut, P., & Heggie, D.C. 2002, MNRAS, 336, 1069 
Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 
Blakeslee, J. P. 1999, AJ, 118, 1506 

Cantiello, M., Blakeslee, J. P., & Raimondo, G. 2007, ApJ, 668, 209 
Cockcroft, R., Harris, W.E., Wehner, E.M.H., Whitmore, B.C., & 

Rothberg, 13. 2009, AJ, paper submitted 
DeCraafT, R.B., Blakeslee, J. P., Meurer, G.R., & Putman, M. 2007, 

ApJ, 671, 1624 

Dirsch, B., Richtler, R., Geisler, D., Forte, J.C., Bassino, L.P., & 

Gicren, W.P. 2003, AJ, 125, 1908 
Evstigneeva, E.A., Drinkwater, M.J., Peng, C.Y., Hilker, M., 

DePropris, R., Jones, J.B., Phillips, S., Gregg, M.D., & Karick, 

A.M. 2007, AJ, 136, 461 
Evstigneeva, E.A., Gregg, M.D., Drinkwater, M.J., & Hilker, M. 

2007 AJ 133 1722 
Fan, Z.', Ma, J., de Grijs, R., & Zhou, X. 2008, MNRAS, 385, 1973 
Fisher, D., Franx, M., & lUingworth, G. 1995, ApJ, 448, 119 
Forte, J.C., Faifer, F., & Geisler, D. 2005, MNRAS, 357, 56 
Forte, J.C., Faifer, F., & Geisler, D. 2007, MNRAS, 382, 1947 
Forte, J.C., Geisler, D., Ostrov, P.G., Piatti, A.E., & Gieren, W. 

2001, AJ, 121, 1992 

Geisler, D., Lee, M.G., & Kim, E. 1996, AJ, 111, 1529 
Georgiev, LY., Goudfrooij, P., Puzia, T.H., & Hilker, M. 2008, AJ, 
135, 1858 

Georgiev, LY., Puzia, T.H., Hilker, M., & Goudfrooij, P. 2009, 

MNRAS, 392, 879 
Gomez, M., & Woodley, L.A. 2007, ApJ, 670, L105 
Harris, W.E. 1996, AJ, 112, 1487 

Harris, W.E. 2001, in Star Clusters, Saas-Fee Advanced Course 28 

(New York: Springer), ed. L. Labhardt & B. Binggeli 
Harris, W.E. & Harris, G.L.H. 2002, AJ, 123, 3108 
Harris, W.E., Harris, G.L.H. , Holland, S.T., & McLaughlin, D.E. 

2002, AJ, 124, 1435 

Harris, W.E., Harris, G.L.H., & McLaughlin, D.E. 1998, AJ, 115, 
1801 

Harris, W.E., Pritchet, C.J., & McClure, R.D. 1995, ApJ, 441, 120 
Harris, W.E., Whitmore, B.C., Karakla, D., Okon, W., Baum, 

W.A., Hanes, D.A., & Kavelaars, J.J. 2006, ApJ, 636, 90 
Harris, W.E., Kavelaars, J. J., Hanes, D.A., Pritchet, C.J, & Baum, 

W.A. 2009, AJ, 137, 3314 
Hasegan, M. et al. 2005, ApJ, 627, 203 
Jordan, A. 2004, ApJ, 613, L117 
Jordan, A. et al. 2005, ApJ, 634, 1002 
King, I. 1962, AJ, 67, 471 

Kobayashi, C, & Arimoto, N. 1999, ApJ, 527, 573 
Kundu, A. 2008, AJ, 136, 1013 



to work, the star formation period needs to stretch over 
at l east 10 Myr and perh aps as long as 20 — 30 Myr 
(see lBailin fc Harris! l2009f ). Contemporary observations 
of dense, massive star-forming regions point encourag- 
ingly in that direction, but are not yet definitive. If the 
average formation time in the most massive GCs is as 
short as ~ 5 Myr, for example, then the self-enrichment 
model would be in serious trouble. 

Many other giant and supergiant ellipticals whose glob- 
ular cluster systems have not yet been studied are well 
within reach of the HST and large ground-based tele- 
scopes. These hold out obvious prospects for increasing 
the GCS sample sizes further by large factors. Major 
strides can be made towards understanding this intrigu- 
ing correlation further. 
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TABLE 1 

Brightest Cluster Galaxies Imaged with the ACS 



Galaxy Cluster Redshift 




E{B-I) (m-M)i R^ff 


B, I Exposures 


or Group (km s~^) 




(Ml < -8.4) (arcmin) 


(sec) 



NGC 1407 


Eridanus 


1627 


-22.35 


0.16 


31.96 


1.17' 


2 X 750, 2 X 340 


NGC 3348 


CfA 69 


2837 


-22.13 


0.17 


33.18 


0.44' 


6 X 1200, 4 X 530 


NGC 3258 


Antlia 


3129 


-21.87 


0.20 


33.23 


0.50' 


4 X 1340, 4 X 570 


NGC 3268 


Antlia 


3084 


-21.96 


0.24 


33.27 


0.60' 


4 X 1340, 4 X 570 


NGC 4696 


CenSO 


2926 


-23.31 


0.23 


33.29 


0.74' 


4 X 1360, 4 X 580 


NGC 7626 


Pegasus I 


3405 


-22.35 


0.17 


33.57 


0.64' 


6 X 1310, 2 X 1300 



TABLE 2 

Characteristics of the Point-Spread Functions 



Galaxy FWHM(S) iV,(S) FWHM(/) Ar,(i") 



NGC 1407 


0.091" 


10 


0.090" 


5 


NGC 3348 


0.081" 


8 


0.089" 


10 


NGC 3258 


0.098" 


6 


0.091" 


18 


NGC 3268 


0.096" 


8 


0.084" 


15 


NGC 4696 


0.102" 


25 


0.095" 


52 


NGC 7626 


0.102" 


6 


0.096" 


11 


Mean 


0.095" 


10 


0.091" 


19 




±o.no.'^3" 




io.oois" 





TABLE 3 

Median Sizes of Globular Clusters versus 
Luminosity 



Mj Range n{blue) (rf^){blue) n{red) {rh){red) 
(pc) (pc) 



(-12.0, -12.7) 
(-11.5, -12.0) 
(-11.0, -11.5) 
(-10.5, -11.0) 
(-10.0, -10.5) 
(-9.5, -10.0) 
(-9.0,-9.5) 

MilkA- W:w 



10 


4.63 


17 


2.21 


105 


3.10 


186 


2.80 


301 


2.65 


419 


2.51 


289 


2.80 



IS 3,06 



15 


3.86 


32 


2.85 


81 


2.63 


150 


2.36 


263 


2.21 


313 


2.21 


139 


2.65 



18 2.10 



TABLE 4 

Median Sizes of Globular Clusters versus Rqc 



Rac/Reff n{blue) {rh}(blue) n{red) {rh)(red) 
(pc) (pc) 



(0.0 - 0.5) 


50 


2.28 


56 


1.95 


(0.5 - 1.0) 


193 


2.51 


185 


2.06 


(1.0 - 1.5) 


232 


2.50 


175 


2.06 


(1.5 - 2.0) 


227 


2.36 


161 


2.18 


(2.0 - 2.5) 


221 


2.63 


141 


2.36 


(2.5 - 3.0) 


159 


2.95 


127 


2.33 


(3.0 - 3.5) 


119 


2.95 


55 


2.36 


(3.5 - 4.0) 


56 


2.95 


32 


2.58 


(4.0 - 5.0) 


37 


3.15 


31 


2.77 
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TABLE 5 

RMIX BiMODAL Fits to the Individual Galaxies 



I range 


Af(bin) 




/"I 












^2 






T2 




Pi 


NOP 1407 






























91 — 99 n 


113 




.885 i 0. 


058 


0. 


.102 it 


.033 


2 


.201 i 0. 


.055 


0. 


120 i 


,033 


0.413 


99 n 9*^ n 






7QQ J- n 


uzo 


n 
u. 


1 oQ -L n 
-Lzo m u. 


ni 7 


2 


225 ^ 


ni ^ 

.U-LO 


n 
u. 


1 91 ^ n 

IZ 1 JZ u 


ni 1 


U.OU 1 


9^4 n — 9^ 


219 




.743 i 0. 


018 


0. 


078 it 


.016 


2, 


177 i 


.020 


0. 


159 i 


,016 


0.298 


9^^ — 9zL n 


216 




.733 i 0. 


.109 


0. 


.109 it 0. 


.018 


2, 


.191 it 0. 


.016 


0. 


158 it 


.015 


0.331 
































21.0 - 23.0 


96 


1 


.864 ±0, 


,061 


0. 


142 it 0, 


.035 


2, 


157 ±0, 


.060 


0. 


104 ±0, 


,031 


0.663 


23.0 - 24.0 


214 


1 


.640 ± 0, 


,012 


0. 


,074 it 0, 


.011 


2, 


.085 ± 0, 


.015 


0. 


,148 ±0, 


.012 


0.318 


24.0 - 25.0 


348 


1 


.642 ± 0. 


Oil 


0. 


,078 it 0, 


.010 


2. 


.085 ± 0, 


.015 


0. 


165 it 0, 


.012 


0.313 


NGC 3258 






























21.0-22.5 


91 


1 


.836 ± 0. 


.032 


0. 


.097 it 0. 


.021 


2. 


.130 ±0. 


.059 


0. 


112 it 0, 


.035 


0.614 


22.5 - 23.0 


113 


1 


.738 ± 0. 


,024 


0. 


.063 it 0. 


.022 


1. 


.997 ± 0, 


.039 


0. 


176 it 0, 


.020 


0.267 


23.0 - 23.5 


197 


1 


.778 ± 0. 


.025 


0. 


.090 it 0. 


.017 


1. 


.133 ±0, 


.037 


0. 


134 it 0, 


.025 


0.483 


23.5 - 24.0 


279 


1 


.753 ± 0. 


.013 


0. 


.085 it 0. 


.009 


2. 


.132 ±0. 


.025 


0. 


137 it 0, 


.019 


0.534 


24.0 - 24.5 


370 


1 


.740 ± 0. 


.009 


0. 


.090 it 0. 


.007 


2. 


.168 it 0, 


.020 


0. 


137 it 0, 


.013 


0.550 


24.5 - 25.0 


402 


1 


.740 ± 0. 


.010 


0. 


.097 it 0, 


.008 


2. 


.146 it 0, 


.015 


0. 


131 it 0, 


.012 


0.513 



NGC 3268 



21.0 


- 22.5 


84 


1 


.815 


±0, 


,017 


0, 


,026 


it 


0, 


,018 


2, 


.050 


±0, 


,023 


0. 


154 ±0, 


.015 


0, 


.128 


22.5 


- 23.0 


101 


1 


.799 


±0, 


,019 


0, 


,030 


± 


0, 


,066 


2, 


,050 


±0, 


,046 


0. 


170 it 0, 


.028 


0, 


.107 


23.0 


- 23.5 


162 


1 


.759 


±0, 


,014 


0, 


,045 


± 


0, 


,013 


2, 


,097 


±0, 


,022 


0. 


198 it 0, 


.015 


0. 


.156 


23.5 


- 24.0 


249 


1 


.768 


±0, 


,010 


0, 


,081 


± 


0, 


,008 


2, 


187 


±0, 


,019 


0. 


152 ±0, 


.015 


0, 


.471 


24.0 


- 24.5 


334 


1 


.878 


±0, 


,014 


0, 


,096 


± 


0, 


,010 


2, 


183 


±0, 


,024 


0. 


165 ± 0, 


.018 


0, 


,449 


24.5 


- 25.0 


377 


1 


.768 


±0, 


.015 


0, 


.108 


it 


0, 


,011 


2, 


.197 


±0, 


,024 


0. 


,163 it 0, 


.017 


0, 


.477 


NGC 4696 








































21.0 


-22.5 


127 


1 


.907 


±0, 


.018 


0. 


.100 


it 


0, 


,000 


2. 


.223 


±0. 


,030 


0. 


150 it 0, 


.000 


0. 


.517 


22.5 


-23.0 


192 


1 


.901 


±0, 


.046 


0. 


.084 


± 


0, 


,027 


2. 


.170 


±0, 


,065 


0. 


,141 it 0, 


.031 


0. 


.378 


23.0 


-23.5 


297 


1 


.846 


±0, 


.010 


0. 


.074 


± 


0, 


,008 


2. 


.220 


±0, 


,018 


0. 


152 it 0, 


.014 


0. 


.389 


23.5 


-24.0 


374 


1 


.814 


±0, 


.009 


0. 


.081 


it 


0, 


,007 


2. 


.234 


±0. 


,015 


0. 


153 it 0, 


.012 


0. 


.441 


24.0 


-24.5 


636 


1 


.812 


±0, 


.010 


0. 


.091 


it 


0, 


,007 


2. 


.213 


±0. 


,019 


0. 


174 it 0, 


.013 


0. 


.445 


24.5 


-25.0 


658 


1 


.847 


±0, 


.012 


0. 


.123 


it 


0, 


,008 


2. 


.278 


±0, 


,023 


0. 


156 it 0, 


.016 


0. 


.589 


25.0 


- 25.5 


516 


1 


.801 


±0, 


.026 


0. 


.122 


it 


0, 


,022 


2. 


.149 


±0, 


,053 


0. 


234 it 0, 


.026 


0. 


.319 


NGC 7626 








































21.0 


- 22.5 


76 


1 


.832 


±0, 


,049 


0, 


100 


± 


0, 


,000 


2, 


116 


±0, 


,030 


0. 


150 it 0, 


.000 


0, 


.194 


22.5 


- 23.5 


218 


1 


.808 


±0, 


,018 


0, 


106 


± 


0, 


,012 


2, 


,203 


±0, 


,020 


0. 


137 ±0, 


.015 


0, 


.437 


23.5 


- 24.0 


228 


1 


.801 


±0, 


,023 


0, 


120 


it 


0, 


,017 


2, 


,219 


±0, 


,025 


0. 


132 ± 0, 


.017 


0, 


,494 


24.0 


- 24.5 


323 


1 


.728 


±0, 


,010 


0, 


,074 


it 


0, 


,007 


2, 


158 


±0, 


,015 


0. 


160 ± 0, 


.012 


0, 


,344 


24.5 


- 25.0 


385 


1 


.752 


±0, 


.020 


0, 


112 


± 


0, 


,014 


2, 


184 


±0, 


,028 


0. 


173 ± 0, 


.019 


0, 


,426 


25.0 


- 25.5 


375 


1 


.741 


±0, 


.024 


0, 


.115 


it 


0, 


,017 


2, 


.182 


±0, 


,037 


0. 


,180 it 0, 


.024 


0, 


.464 
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TABLE 6 

RMIX BiMODAL Fits to the Combined Sample 



Mj range A'^(bin) 



M2 



<T2 



Pi 



285 


1 


757 


±0 


038 





102 


±0 


032 


1 


944 


± 





029 





209 


±0 


013 





181 


±0 


109 


359 


1 


642 


±0 


017 





100 


±0 


000 


1 
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Oil 



-12.7, -11.5) 
-11.5,-11.0) 
-11.0, -10.5) 
-10.5,-10.0) 
-10.0, -9.5) 
-9.5, -9.0) 
-9.0, -8.5) 
-8.5, -8.0) 

-12.4. -10.4) 
-10.4,-8.4) 



TABLE 7 
Radial Metallicity Gradients 



Galaxy 


Ar(blue) 


a 


b 


W(red) 


a 


6 


NGC1407 


115 


1.630 ±0.010 


-0.033 ± 0.035 


225 


2.036 ± 0.009 


-0.027 ± 0.028 


NGC3268 


450 


1.597 ±0.007 


-0.042 ± 0.019 


292 


1.980 ± 0.009 


+0.002 ± 0.025 


NGC3258 


386 


1.604 ± 0.008 


-0.033 ± 0.025 


306 


1.970 ± 0.007 


-0.030 ± 0.024 


NGC3348 


141 


1.581 ± 0.012 


-0.017 ±0.035 


148 


1.989 ± 0.011 


-0.036 ±0.038 


NGC4696 


590 


1.628 ± 0.004 


-0.038 ±0.014 


576 


2.010 ± 0.006 


-0.011 ±0.018 


NGC7626 


310 


1.615 ± 0.007 


-0.060 ±0.023 


326 


2.035 ± 0.007 


-0.068 ±0.026 


All 6 


1992 


1.615 ±0.003 


-0.0379 ± 0.0068 


1873 


2.008 ± 0.003 


-0.0371 ± 0.0079 
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Fig. 1. — Upper panel: Distribution of half-light radii for Milky Way globular clusters. All 102 such clusters with luminosities brighter 
than Mv = —6.5 are included here, to roughly match the lum inosity range of the GCs observed in the six target galaxies of this study. 
Ninety percent of them have r^i < 6 parsecs. Data are from lHarrisI 1199^) . Lower panel: Distribution of King-model central concentration 
c = log(rt/rc) for the same sample of Milky Way clusters above. The bin at c = 2.5 is populated by the core-collapsed GCs, for which 
c = 2.5 is arbitrarily assigned. 
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Fig. 2. — Selection of candidate stars for defining the stellar point spread function. The data shown here are SExtractor measurements 
of the half-light radius ro.s versus an aperture magnitude, for all detected objects in the NGC 4696 galaxy field. No attempt has been 
made to cull out nonstellar objects or detections due bad pixels, clumps of noise, or other artifacts. The magnitude scales have arbitrary 
zeropoints and do not represent the true B, I magnitudes. Note in each panel the narrow sequence of "starlike" sources; at the very bright 
end (m/ < 12, nig < 12.5) these become saturated and the sequences curve upward. The dashed boxes mark out the objects that were 
chosen as candidate PSF stars for further investigation (see text). 
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Fig. 3. — Three sample illustrations of the ISHAPE profile fits to objects in the NGC 4696 field. In each panel the upper right quadrant 
is the real object, the upper left quadrant shows the best-fitting model (that is, the KingSO profile convolved with the PSF), and the lower 
left quadrant shows the residuals (real — model). The lower right quadrant shows the pixel weight map. (a) Left panel: A bright nonstcUar 
object: the best-fit solution is for FWHM a = 1.3 px, cUipticity b/a = 0.91, and S/N = 440. In this case the intrinsic size of the object is 
nearly as large as the PSF width of ~ 1.8 px. (b) Middle panel: A faint nonstcUar object, for which a = 0.82 px, b/a = 0.50, and S/N = 24. 
Note that the object looks much rounder than its intrinsic shape {b/a) = 0.5 because it has been convolved with the circularly symmetric 
PSF, whose FWHM is more than twice as wide, (c) A starlike object, for which a = 0, b/a = 1, and S/N = 108. 
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Fig. 4. — Curves of growth for simulated globular clusters. For each of the two graphs, the uppermost curve shows the magnitude 
enclosed within an aperture of radius r for a starlike object. All the other curves in each set show the curve of growth for simulated GCs 
of various sizes, convolved with the same stellar PSF profile. All the curves are plotted to converge at r = 20 px (1 arcsecond). The more 
extended the GC profile is, the more gradually the curve of magnitude versus r increases: from top to bottom, the assumed GC profile 
FWHM values are (0.0,0.2,0.4,0.6,0.8,1.0,1.25,1.50,1.75,2.0,2.25,2.50,2.75,3.0). The vertical bar plotted at r = 2.5 px is the adopted 
radius for our fixed-aperture photometry, and it covers the typical range of GC profile widths for the GCs in our data. 
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Fig. 5. — Aperture eorrections plotted as a function of intrinsic size of the individual globular cluster. These are generated from the 
family of curves plotted in the previous figure, where the magnitude enclosed within r = 20 px is defined as the "total" magnitude of the 
object. To obtain this total magnitude, the magnitude within a 2.5-px aperture is measured and then the aperture correction from this 
graph is subtracted. The data for our actual GCs fall in the regime FWHM < 0.6 px, i.e. the lower left part of the curves. 
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Fig. 6. — Scatter plot of the measured ISHAPE sigiial-to-noisc in B versus the S/N for the same objects in /. Ideal 1:1 agreement would 
fall along the dashed line, while the solid line shows the actual correlation where S/N(B) ~ 0.4 S/N(7). Although the data for only the 
NGC 4696 field are shown, these results are typical of all six of the galaxy fields. 
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Fig. 7. — Tests of the sensitivity of our GC size measurement to the assumed central concentration c of the King model. The panels on 
the left show tests in the NGC 4696 field of the differences in measured rj, relative to the baseline King.30 model, for (Kingl5-King30) and 
(Kingl00-King30). On the right, similar results are shown from the NGC 3268 field for the comparison (KinglO-KingSO) and (KinglOO- 
KingSO) . The units of both axes are in ajrcseconds. 
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Fig. 8. Tests of the sensitivity of our sohitions for GC size to the assumed eUipticity e. The simulated clusters have an assumed KingSO 
profile with FWHM= a, coii\'olvecl witli a model PSF with FWHM=2.0 pixels. In eaeh graph the correct (input) value of the intrinsice 
FWHM of the clusters is shown as the horizontal dashed line. In the left-hand set of graphs, ISHAPE was allowed to solve for both the 
FWHM and the eUipticity, while in the right-hand set, it was forced to assume e = and solve only for the FWHM. 
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Fig. 9. — Tests of ISHAPE solutions on simulated clusters for the object shape ib/a) = (1 — e) and orientation angle Q as functions of 
object brightness and size. The correct (input) values are {b/a) = 0.8 and 9 = 45°, as shown by the horizontal dashed lines. 
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Fig. 10. — Tests of the ishape measurements for simulated globular clusters of various intrinsic sizes. KING30 model GCs of an assumed 
input FWHM (measured in pixel units, where 1 px = 0.05") and a range of luminosities are convolved with the stellar point-spread 
function and inserted into the image on top of realistic background noise. These simulated objects are then measured with ISHAPE and 
the measured FWHM is plotted against tiie input value. The dashed lines show the FWHM of the assumed stellar PSF at 2.0 px. These 
tests show that accurate values of the GC sizes can be recovered down to FWHM ~ 0.2 px, or 10% of the stellar FWHM. 
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Fig. 11. — Comparisons of the ishape measurements of object size for the real objects in each of the six galaxy fields. Here the best-fit 
half-light radius r^i = 1.48 fwhm from ISHAPE is shown in B versus / for all objects with ISHAPE total S/N > 50. Ideal one-to-one 
agreement is shown by the solid lines. 
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Fig. 12. — The color-magnitude data for the measured objects in each of six target galaxy fields. Each individual object is plotted 
with its fully size-corrected {B, I) magnitudes as described in the text. The majority of the objects in each field are the globulcir cluster 
populations. The main source of field contamination is predominantly from small, faint background galaxies. 
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Fig. 13. — The color-magnitude data for all six galaxies, now converted to luminosity and intrinsic color {Mj,{B — I)o) and with all 
fields combined. The left panel shows the size-corrected aperture photomet ry from the current discussion, while the right panel shows the 
data for the same fields based solely on PSF fitting as published earlier in [Hams et al.l pOOeh . The mean lines through the blue and red 
globular cluster sequences for the newer data at left are shown by the heavy dashed lines. The same lines are repeated in the right panel 
to indicate the small offsets between the two datasets. Sample crrorbars showing the measurement uncertainties of the photometry are 
shown for Mj = -8, -9, -10. 
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Fig. 14. — Measured half-light radii for the individual globular cluster candidates in the six galaxy fields. Only the objects with r/j > 1 pc 
are plotted, though they have not been selected by {S/N). The two metallicity groups (blue and red clusters) can be seen in each sample. 
Over all six fields, the total is 6476 objects. 
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Fig. 15. — Histogram (solid broken curve) of the best-quality size measurements {S/N > 50) for 3330 objects as determined by ISHAPE, 
over all six galaxy fields combined. All measurements of the object size are converted to half-light radius in parsecs. The shaded area 
covering rj, < 2 pc indicates the approximate size range within which ISHAPE cannot easily distinguish starlike from nonstellar objects. 
For comparison, the block histogram on the left side of the graph peaking at = shows the distribution of measured radii that ISHAPE 
returned for 5300 simulated starlike objects at an assumed distance of 40 Mpc; all of these by definition have true radii rj, = 0. Finally, the 
dotted histogram shows the size distribution for the 66 globular clusters in the Milky Way that are more luminous than the GC luminosity 
function turnover point and no further than 20 kpc from the Galactic center (see text). 
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Fig. 16. — Left panel: Measured half-light radii for the combined globular cluster populations in all six galaxies. The best-quality 2495 
objects with sizes larger than = 1 pc and with (S/N) > 50 are shown. Right panel: Histograms of the values for the metal-poor GC 
population (1.3 < (B ■— I)o < 1.8, solid line) and the metal-rich population (1.8 < (B — /)o < 2.4 dashed line). 
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Fig. 17. — Upper panel: Half-light radii for the metal-poor globular clusters, plotted versus cluster luminosity M/. Data are the same as 
in the previous figure. Median values of from Table 3 in half-magnitude intervals are shown as the large dots. Lower panel: The same 
data for the metal-rich globular clusters. On average the metal-poor clusters have median sizes about 15% larger than do the metal-rich 
ones, at any luminosity. The heavy dashed lines at the right in each panel show the mean size/luminosity relation for UCDs and dE nuclei, 
from lEvstigneeva et al.l II2008I ). The dotted lines show the lower envelope of the same UCD/dE,N relation. 
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Fig. 18. — The half-light radii for globular clusters plotted as a function of color (metallicity ) , for three radial zones of galactocentric 
distance. 
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Fig. 19. — The half-light radii for globular clusters plotted as a function of galactocentric distance. Upper panel: Blue, metal-poor GCs 
defined as in the previous figures; Lower panel: Red, metal-rich GCs. The large connected dots show the empirical trend of the median 
value. The dashed lines give the power-law fits to the median points, with equations as given in the text. 
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Fig. 20. — Deconvolution of the (B — I) color distributions for the GC populations in the six individual galaxies. The fitting code RMIX 
(see text) has been used to isolate the best- fit mean colors of the blue and red sequences, as listed in Table [5] 




Fig. 21. — Mean points in (M/, {B — I)o) for all six galaxies, as determined from the bimodal Gaussian fits listed in Table[5] Solid dots 
are the fiducial points for the blue sequence, open triangles for the red sequence. The solid lines through each sequence connect the mean 
points as determined from the combined dataset in Table |6] 
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Fig. 22. — Left panel: CMD for all objects with positive measurements of effective radius Vf^. Right panel: CMD for all objects with 
either "starlike" size measurements (r^i = 0) or those too faint for successful ISHAPE fits. The heavy dashed lines showing the red and 
blue GC sequences are the same as in previous figures. 
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Fig. 23. — Le/t panel: CMD for all stailikc objects in the outer regions of the three "control field" galaxies, NGC 5322, 5557, and 
7049. The data have been measured and converted to absolute magnitude and color as described in the text. Right panel: CMD for all 
control-field objects with positive size measurements. Note the presence of globular clusters along the normal blue and red sequences in 
the outskirts of these three galaxies. 
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Fig. 24. — Histograms in (B — I)o for the measured objects in all the galaxy fields combined, subdivided by GC luminosity Mj. The 
unshaded histograms show the numbers of objects per 0.05-mag color interval in the six program galaxies. The shaded histograms show the 
numbers in the composite "control field" (see text), multiplied by a factor of 6 to match the total areas. The control field data represent 
only an upper limit to the true background contamination level. 
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Fig. 25. — Upper panel: color-magnitude distribution for the GCs in all six galaxies, in the luminosity interval Mj < —10.4. Lower 
panel: color-magnitude distribution for the fainter objects, Mj > —10.4. The reader is invited to block off each panel in turn to view the 
other one independently. 
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Fig. 26. — Deconvolution of the (B — I) color distributions for the GC populations in the combined data sample. The upper panel shows 
the histogram and the RMIX two-gaussian fit for the brightest clusters, in the luminosity range Mj = (—12.4, —10.4). The lower panel 
shows the same data for the less luminous clusters in the range Mj = ( — 10.4, —8.4). Note that the faint limit Af/ = —8.4 is approximately 
the level of the GCLF turnover point (see text). 
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Fig. 28. — CMD for the 7831 measured objects with size measurements > 1.5 pc, drawn from the previous Fig. 1221 This sample 
minimizes field contamination and provides the best visual representation of the blue and red GC sequences. The heavy dashed lines show 
the mean points along the two sequences, as listed in Table |6] The solid lines shows the polynomial fits to these mean points as given 
in equations (5,6), and the dotted lines show the same polynomial fits determined directly from the individual cluster data divided at 
{B - I)o = 1.8 (see text). 
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Fig. 29. — Left panel: Combined color-magnitude data for tlie six BCGs, from Fig. 13a. The histogram at top shows the distribution in 
(B — /)o for the luminosity range —10.0 < Mj < —8.5, a range fainter than the main region of the MMR but brighter than the region where 
the colors are broadened significantly by observational measurement scatter. Right panel: Simulated data for a bimodal GC population 
with parameters as described in the text. Both the red and blue sequences are purely vertical (no MMR), with mean colors {B — I)o = 1-56 
and 2.00. The histogram above shows the color distribution for the same luminosity region (—10.0, —8.5). 
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Fig. 30. — Left panel: Color-magnitude distribution for the GCs with S/N > 20 and measured half-light radii smaller than 2 pc. Middle 
panel: The same data for objects between 2 and 4 pc. Right panel: The same data for the biggest objects, > 4 pc. In all three panels, 
the dashed lines show the mean GC sequences derived from the entire data sample (see Fig. 1131 and accompanying text). 
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Fig. 31. — Color-magnitude distribution for the globular clusters with S/N > 20 in all six galaxies, subdivided by radial zone. The 
projected galactocentric distance of each object is expressed in units of the effective radius Reff of the galaxy's spheroid light. In each 
panel, the dashed lines show the mean GC sequences derived from the entire data sample. The most remote measured objects are at radii 
~ 5Rsff (see text). 
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Fig. 32. — Mean points determined from the RMIX bimodal-Gaussian fits for the blue cluster sequence in three radial zones as defined 
in the previous figure. Solid dots and line are for the inner zone, the triangles and dashed line for the middle zone, and the crosses and 
dotted line for the outer zone. 
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Fig. 33. — Fraction of luminous, blue-sequence globular clusters bluer than (B — 7)o = 1.8, plotted against color. Dotted line is for the 
inner-zone clusters {R < R^ff), dashed line is for the mid-zone {Reff < R < ^Reff)^ a^id solid line is for the outer zone {R > 2Rgff). 



56 



Harris 




Fig. 34. — Colors of the bright (M/ < —9) objects with positive size measurements {rf^ > 0), plotted against galactocentric distance. 
Data from all six galaxies are combined. The individual distances R are normalized to the effective radius R^jtoi the spheroid light. The 
solid lines show the best-fit power-law solutions for the mctallicity gradients, as listed in the last line of Table |7] 



